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Abstract 


Microporous polyvinylidene fluoride (PVDF) membranes were prepared by casting and coag- 
ulating solutions of the polymer in eight different solvents. The membranes were characterized 
by scanning electron microscope observations, water content measurements and water permea- 
bility tests. The possible correlations between the resulting membrane characteristics and the 
properties of both the polymer-solvent solutions and the solvent-nonsolvent systems are exten- 
sively discussed. The analysis indicates that most of the correlations suggested or adopted in the 
literature do not hold for PVDF membranes, whereas the transport properties of the solvent- 
nonsolvent system seem to affect significantly the final membrane structure and the resulting 
performances in separation processes. 


Keywords: membrane preparation and structure; microporous membrane; polyvinylidene fluoride 
membrane; scanning electro microscopy; water content 





Introduction 


Phase separation is a widely used method for making microporous polymer 
membranes. In this process a viscous solution is first cast on a suitable support, 
then immersed in a nonsolvent bath. The membrane is formed by polymer 
precipitation, which occurs as a consequence of concentration variations fol- 
lowing diffusive interchange between the solvent and the nonsolvent. The rate 
of polymer precipitation is determined at each point by the progress of the 
concentration, which is in turn governed by the interchange rate. 

Depending on the rate of polymer precipitation, the following three types of 
membranes can be obtained: 

(a) symmetric, with an almost even porosity across the membrane cross- 
section; 

(b) asymmetric, with a selective thin microporous upper layer (skin) on a 
thicker macroporous globular or spongy sub-layer; 
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(c) asymmetric, with large voids and/or finger-like cavities beneath the mi- 
croporous upper layer. 

There is evidence that a low precipitation rate leads to type (a) membranes, 
whereas a high precipitation rate gives type (c) membranes. Nonetheless 
membrane formation is complex to describe, as it occurs in very short times 
and involves a great number of elementary processes. In consequence, the mo- 
delling of membrane formation is very difficult for three main reasons: (i) the 
number of the intervening elementary processes is very high; (ii) many in- 
volved mechanisms and phenomena (such as convective contributions, ther- 
mal effects, demixing kinetics, dependence of the chemical potential on the 
system composition and the stress distribution) are still poorly described or 
understood; (iii) the independent evaluation of many relevant parameters (in- 
cluding ternary diffusivities, interaction parameters, and boundary and initial 
conditions) is difficult or highly approximate. However, it is likely that in some 
systems a few or just one mechanism could control the process. In such cases 
only some or one parameter could be important, so that it could be worthwhile 
to look for some consistent correlation between a single parameter and the 
obtainable membrane structure. This alternative approach can furnish both 
practical guidance in membrane preparation and an insight into the process. 
Of course this approach is not intended to substitute the modelling of the pro- 
cess; preferably, the holding correlation should be related back to a possible 
model and justified by it. 

In this sense, several experimental observations and qualitative comments 
on the effect of polymer-solvent and/or solvent—nonsolvent interactions on 
the structure and properties of the resulting membrane are reported in the 
literature [1-17]. A large amount of work has been carried out on a limited 
class of polymers, essentially cellulose acetate and its derivatives [1-—5,7,9- 
11,14-16], polyamides [4,13,15,17] and polysulphones [10,12,14,15], but very 
little information is available on polyvinylidene fluoride (PVDF), which now- 
adays is widely used as a membrane material since it exhibits excellent chem- 
ical resistance and good physical and thermal stability [18]. 

The purpose of this work is to investigate for the PVDF system the existence 
of correlations between the properties of both PVDF-solvent solutions and 
solvent—nonsolvent systems and the structure and properties (water content, 
water flux, thickness and void volume) of the membranes. An attempt is also 
made to draw general conclusions and inferences, so that the obtained indi- 
cations should furnish valuable help in the analysis of other systems and in the 
development and improvement of models of the process. 


Experimental 


Materials 
The PVDF was a commercial product (Foraflon 1000 HD, Ugine Kuhl- 
mann) with an average molecular mass M=4.5x 10°. The total polymer sol- 











TABLE 1 


Density, viscosity, solubility parameters and mutual water diffusivity of PVDF solvents 





Solvent ps" ns* Jas Ops Ons d,s” Ds _wX 10° Dw 5X 10°D,,, x 10° 
(kg/m*) (mPa-sec) (MPa'/*) (MPa'’?) (MPa'’?) (MPa'’”) (cm?/sec) (cm/sec) (cm*/sec ) 








DMA 941.2 0.9472 16.8 11.5 10.2 22.7 9.1 16.8 11.8 
DMF 949.1 0.8499 17.4 13.7 11.3 24.8 10.2 17.1 12.8 
DMSO 1100.4 2.1878 18.4 16.4 10.2 26.7 10.7 6.9 8.4 
HMPA 1025.8 3.5570 18.4 8.6 11.3 23.2 6.2 6.4 6.3 
NMP 1032.4 1.8179 18.0 12.3 UP : 22.9 8.9 9.3 9.1 
TEP 1069.4 1.6753 16.8 11.5 9.2 22.3 6.3 13.7 8.7 
TMP 1213.4 2.1937 16.8 16.0 10.2 22.3 8.0 9.2 8.6 
TMU 968.1 = 1.5330 16.8 8.2 11.1 21.7 7.8 12.0 9.5 











“Temperature 20°C; Ref. [20]. 
»Temperature 25°C; Ref. [21], pp. 153-158. 


ubility parameter was 6, p= 23.2 MPa'’’ and the dispersion (d,,p), polar (0, p) 
and hydrogen bonding (0,,») parameters were 0g p~=17.2, 0,p=12.5 and 
Oy p=9.2 MPa’’’, respectively [19]. 

The polymer solvents were of reagent grade: N,N-dimethylacetamide 
(DMA), N,N-dimethylformamide (DMF), dimethylsulphoxide (DMSO), 
hexamethylphosphoramide (HMPA), N-methyl-2-pyrrolidone (NMP), te- 
tramethylurea (TMU), triethyl phosphate (TEP) and trimethyl phosphate 
(TMP). Some physical constants such as density (pg), absolute viscosity (5 ) 
and solubility parameters for these solvents are listed in Table 1. 


Properties of PVDF -solvent solutions 

Limiting viscosity numbers at 20°C of PVDF in each solvent were taken 
from Ref. [19]. The kinematic viscosity (vcs) of 15 wt.% casting solutions was 
measured at 20°C by using a Cannon Fenske capillary viscometer. Absolute 
viscosity (cs) is defined as: 


Ncs =PcsY’cs (1) 


where the solution density pes was measured with a pycnometer. Specific vis- 
cosity (4% ) was determined by the ratio: 

nun «m Shin 
nes = ‘ics = Ns (2) 


Ns 


Properties of water-solvent systems 
Excess enthalpies of mixing at 37°C of the water-solvent systems were taken 
from Ref. [20]. 


Mutual diffusivities Ds yw of solvent at very low concentrations in water and 
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Dw_s of water at very low concentrations in solvent were evaluated from the 

Wilke—Chang equation [22]: 

(~,M,)'/°T 
mVa° 





D,.»=7.4X107" (3) 
where D,, ;, (cm”/sec) is the liquid mutual diffusivity of a in essentially pure b; 
T (K) is the absolute temperature; M is the molecular mass; V (cm*/mol) is 
the molar volume; 7 (Pa-sec) is the absolute viscosity; and 9 is the association 
factor. 

The following assumptions are made for computations: 
(i) the molar volume of water is taken as 4 times greater than the actual value 
as suggested in [23]; 
(ii) association factors are 2.26 and 1.1 for Ds w and Dw sx, respectively. 


Cloud points 

A set of ternary system samples with increasing water concentration was 
prepared by adding known amounts of water-solvent mixtures of fixed com- 
position to weighed PVDF solutions. The samples were continuously stirred at 
50°C for 2 hr, then at 20°C for 8 hr. The cloud point was assumed to be at the 
concentration of the first sample of the set showing constant turbidity at visual 
inspection. The solubility diagram was obtained by repeating this procedure 
on several sets with different solvent contents. 


Membrane preparation 

Solutions of different PVDF concentration were cast at 20°C ona glass plate 
as films of ~ 350 um thickness, and after 30 sec, immersed in water at 4°C. 
Membranes were leached for 48 hr under running water before characterization. 


Membrane characterization 

Membrane cross-sections were prepared according to the procedure reported 
in Ref. [24] and observed by a Cambridge Stereoscan 250 MK2 scanning elec- 
tron microscope (SEM). 

Membrane water content was determined as follows: 


w- fe) 


(4) 


where W and W, are the weights of the wet and dried membrane, respectively. 

The degree of crystallinity of the polymer in the membrane was measured 
by a calorimetric method in order successively to evaluate the polymer density 
in the membrane, as the densities of amorphous and crystalline PVDF are 1680 
and 1920 kg/m*, respectively [19]. A differential scanning calorimeter (Per- 
kin-Elmer DSC IIB) was used, and the enthalpy of fusion of PVDF crystals 
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was assumed to be 6.7 kJ/mol, according to Ref. [25]. The degree of crystal- 
linity ranged from 55-60% according to the type of solvent used, and conse- 
quently the polymer density obtained in the membrane was approximately 
1820 kg/cm’®. 

The knowledge of the polymer density in the membrane and of the value of 
the water content allowed assessment of the void fraction via the hypothesis 
that in the wet membrane all the voids are filled with water. 

The thickness of the wet membrane was measured with a micrometer. The 
product of the thickness of the wet membrane and the void fraction gives a 
measure of the volume of voids or cavities per unit surface area. 

Transmembrane water flux was evaluated with a 2 hr permeation run at 
20°C for a pressure difference of 200 kPa and a 5 m/sec recirculation rate. 
Details of the permeation unit are reported elsewhere [26]. 


Results and discussion 


Membrane structure and properties 

Figure 1 shows SEM micrographs of cross-sections of membranes cast from 
15 wt.% PVDF solutions. Most of the membranes exhibit highly inhomoge- 
neous structures owing to the presence of large voids and cavities of different 
size and shape beneath the upper layer. SEM micrographs reported in Figs. 2- 
9 supply more detailed information on the structure of the membrane through 
magnified views of the upper and the lower parts of the membrane cross-sec- 
tion. By using TEP as solvent a membrane forms with a honeycomb structure 
without cavities (Fig. 2). The size of the alveoli progressively increases from 
the top to the bottom of the membrane. The DMF cast membrane (Fig. 3) 
exhibits a large number of short finger-like cavities beneath the upper layer. 
The cavities are separated by walls composed of discrete polymer globules. The 
structure of the porous sublayer is globular and the size of globules appears to 
be fairly uniform over the entire membrane cross-section. In the membranes 
cast from DMA (Fig. 4), TMU (Fig. 5) and TMP (Fig. 6) the cavities pro- 
gressively extend through a considerable fraction of the membrane thickness. 
In addition, it is apparent that the globules near the upper layer are more 
compacted and larger than those near the bottom. By using NMP (Fig. 7) or 
DMSO (Fig. 8) as solvent, cavities are seen to grow in breadth and length 
towards the bottom of the membrane without a well defined shape. In the mem- 
brane cast from HMPA (Fig. 9) short cavities and isolated voids appear in the 
porous globular sublayer. 

An increase of the polymer concentration in the casting solution leads to a 
reduction of both the number and the size of cavities, but the structure always 
remains of the globular type. This effect is apparent in Figs. 10 and 11, where 
SEM micrographs of cross-sections of membranes cast from NMP and TMP 
solutions, respectively, of different PVDF concentrations are reported. 











TMP 


TEP 


NMP 


Fig. 1. Scanning electron micrographs of cross-sections of membranes cast from 15 wt.% PVDF 


solutions in different solvents. 
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Fig. 3. Scanning electron micrographs of cross-sections of the membrane cast from 15 wt.% PVDF 
solution in DMF: (a) overall cross-section; (b) upper part of (a); (c) bottom part of (a). 








Fig. 4. Scanning electron micrographs of cross-sections of the membrane cast from 15 wt.% PVDF 
solution in DMA: (a) overall cross-section; (b) upper part of (a); (c) bottom part of (a). 


Membrane properties are listed in Table 2. 

Of all the investigated solvents HMPA yields membranes with the highest 
water content and water flux. It is worthwhile to observe that an increase in 
both the number and the size of the cavities in the sublayer does not necessarily 
cause an increase in the water flux. This is because the water flux is mainly 
controlled by the porosity of the upper layer. Indeed, preliminary results ob- 
tained from porosity measurements and ultrafiltration tests with aqueous so- 
lutions of dextran (nominal average molecular mass M ~ 10,000; Fluka Chemie 
AG) indicate that the membrane cast from HMPA (the structure of which 
appears more homogeneous than that of the membranes cast from NMP or 














Fig. 5. Scanning electron micrographs of cross-sections of the membrane cast from 15 wt.% PVDF 
solution in TMU: (a) overall cross-section; (b) upper part of (a); (c) bottom part of (a). 








Fig. 6. Scanning electron micrographs of cross-sections of the membrane cast from 15 wt.% PVDF 
solution in TMP: (a) overall cross-section; (b) upper part of (a); (c) bottom part of (a). 


DMSO) shows a higher number of larger pores in the upper layer. In conse- 
quence it exhibits a greater permeate flux and a much lower dextran retention. 

When the different membranes are sorted according to any property, the 
resulting sequence is always the same. This is shown in Fig. 12, where nor- 
malized values of the properties are plotted for the different membranes. This 
fact agrees with the consideration that the morphology of the sublayer in an 
asymmetric membrane (which is somehow described by the values for the water 
content, membrane thickness and cavity volume) should be mainly deter- 
mined by the structure of the skin [3,13], which the value of the water depends 
on directly. Actually, the skin should be formed in very short times (very small 
fractions of a second) and consequently acts as a barrier against subsequent 











Fig. 7. Scanning electron micrographs of cross-sections of the membrane cast from 15 wt.% PVDF 
solution in NMP: (a) overall cross-section; (b) upper part of (a); (c) bottom part of (a). 





Fig. 8. Scanning electron micrographs of cross-sections of the membrane cast from 15 wt.% PVDF 
solution in DMSO: (a) overall cross-section; (b) upper part of (a); (c) bottom part of (a). 


interdiffusion of solvent and nonsolvent in the bulk of the forming membrane. 
Thus the structures of the skin and the underlying zone are correlated to some 
extent for the investigated asymmetric membranes. As a consequence, for PVDF 
membranes the water flux should not be considered merely as a measure of the 
membrane permeability, which is directly determined by the skin porosity, but 
can be used also as an indirect indication of sublayer morphology. 


Factors affecting membrane structure and properties 

Many authors have searched for correlations between membrane structure 
and one or some properties of the system used to form the membrane. For this 
reason, the properties of the current eight water-solvent-PVDF systems have 
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Fig. 9. Scanning electron micrographs of cross-sections of the membrane cast from 15 wt.% PVDF 
solution in HMPA: (a) overall cross-section; (b) upper part of (a); (c) bottom part of (a). 





Fig. 10. Scanning electron micrographs of cross-sections of membranes cast from 10 wt.% (a), 15 
wt.% (b) and 20 wt.% (c) PVDF solutions in NMP. 


been investigated in order to determine the existence of correlations, to check 
the correlations suggested in the literature, and simultaneously to give some 
insight into or knowledge of the mechanism of membrane formation. 


(i) Polymer-solvent interactions 

The effect of polymer-solvent interactions on the structure and properties 
of the membrane has been studied by various authors [2,7,8,10,11,14]. It is 
generally observed that the lower are these interactions, the higher is the rate 
of polymer precipitation, with consequent formation of finger-like structures. 
Strathmann [8,13,14] assumes as a measure of the compatibility of polymer 
and solvent the difference of their solubility parameters: 


Ov.p-s=O.p — Os (5) 











TABLE 2 


Fig. 11. Scanning electron micrographs of cross-sections of membranes cast from 10 wt.% (a), 15 
wt.% (b) and 20 wt.% (c) PVDF solutions in TMP. 





Properties of membranes cast from 15 wt.% PVDF solutions in different solvents 








Solvent Water flux Water content Thickness of wet Void volume per 
(m/sec ) (%) membrane (sim) unit area (um) 
DMF 0.7 52 135 88 
DMA 1.7 57 150 104 
TMU 0.5 60 150 108 
TEP 2.7 60 140 101 
NMP 11.6 67 200 155 
TMP 11.6 67 220 171 
DMSO 24.3 67 240 186 
HMPA 38.2 81 270 237 





and concludes that the larger this difference is, the higher is the tendency to 
form finger-like structures. By comparing the values of 0, p_s for the PVDF- 
solvent systems reported in column 2 of Table 3 with the SEM micrographs 
shown in Figs. 1-9, it is possible to conclude that the membrane structure does 
not depend on 6, p_s. Due to the particular chemical nature of PVDF and the 
solvents, it is preferable to consider a more sophisticated formulation where 
the three components dg, 6,, and 6), of 6, appear. According to Hansen [27]: 


dj =d4 +0; + di, 


and hence 


d.,p-s = [ (dap —9a,s)* + (bp.p —5p,5)? + (One — Ons)? ] 


Here 0, p_s represents the distance between the polymer and the solvent coor- 


(6) 


(7) 
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Fig. 12. Bar graph of four normalized membrane properties as a function of the solvent used. 


TABLE 3 


Solubility parameter disparity of PVDF and solvents and limiting viscosity numbers of the polymer 








Solvent O:.p-_s* b..p-s° [7] 
(mPa’/*) (mPa’’*) (ml/g) 
DMA —0.5 1.5 129 
DMF 1.6 2.4 117 
DMSO 3.5 4.2 105 
HMPA 0 4.6 158 
NMP —0.3 2.2 128 
TEP —0.9 1.1 130 
TMP —0.9 3.7 115 
TMU —1.5 4.7 142 





“From egn. (5). 
>From eqn. (7). 


dinates in a three-dimensional dq, 6,, 6, Hansen space. The lower this distance, 
the better the solvent. The 6, p_s values calculated from eqn. (7) are reported 
in column 3 of Table 3. But even this approach yields no significant correlation 
between the values of the solubility parameters of the casting solution com- 
ponents and the structure of the membrane. 

Strictly speaking, the quality of the solvent is determined by the value of the 
osmotic second virial coefficient, A, [28], since the higher the value of this 
coefficient the better the solvent. For solutions of flexible-chain polymers, such 
as PVDF [29], a correlation holds between A, and [7]. By comparing the val- 
ues of [7] for the PVDF-solvent systems reported in the last column of Table 
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3 with the SEM micrographs in Figs. 1-9, one can conclude that no correlation 
exists between polymer-solvent compatibility and the membrane structure. 

For concentrated solutions, Frommer et al. [2] indicate that the concentra- 
tion of nonsolvent (water) at which the polymer precipitates is somehow a 
measure of the solvation power, which can be correlated with the membrane 
structure. Furthermore, this quantity can furnish an appraisal of solvent—non- 
solvent interactions. Cloud point data for PVDF-solvent—water systems are 
reported in the ternary phase diagram shown in Fig. 13. The amount of water 
required for precipitating PVDF increases, at a fixed PVDF concentration, 
with the following order for the different solvents: 


HMPA>TMU>NMP>DMA>DMF>TEP>TMP>DMSO 


However, even in this way the order does not correspond to that found for the 
membrane structures obtained from the different solvents. 


(it) Solvent-nonsolvent interactions 

It has been shown by various authors [3,5,8-13] that a high tendency of a 
solvent to mix with the nonsolvent plays an important role in membrane po- 
rosity and the formation of finger-like cavities. So et al. [4] report on the 
preparation of asymmetric membranes cast from several polymer-solvent so- 
lutions immersed in water, and find, in most cases, that the higher the differ- 
ence of the solubility parameter of water (0, yw=47.8 MPa'’* from Ref. [21] ) 
and solvent 


~ 


Ouw s =0.w — drs (8) 


and hence the lower the tendency of mixing [6], the higher the membrane 
water content. Values of 0, w_s for the present systems can easily be obtained 


DMSO 
TMP 
TEP 
OMF 
DMA 
NMP 
9 v TMU 
| 9 HMPA 
v 


®oo 8 »e 


D> 








20 30 H20 


Fig. 13. Ternary phase diagram at 20°C for PVDF -solvent-water systems. 
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once 0, s is read from Table 1. However, the calculated values appear to be 
random with respect to the water content of the PVDF membranes in Table 2. 
Frommer and Lancet [3] assume the enthalpy of mixing to be an indicative 
parameter of the tendency of water and of the solvent to dissolve mutually in 
one another, and show that high values of this parameter correspond to high 
rates of polymer precipitation, which favour the formation of finger-like mem- 
brane structures. Similar trends are reported in Refs. [5], [8] and [13]. The 
values of the excess enthalpy of mixing for the investigated water-solvent sys- 
tems [20] are plotted as a function of composition in Fig. 14. Again, by com- 
paring this figure with Figs. 1-9, it can be concluded that on the basis of the 
thermal effect alone it is not possible unequivocally to predict the final struc- 
ture of the PVDF membranes. Actually, the tendency of mixing is more pre- 
cisely described by the free energy of mixing. However, the influence of this 
quantity has already been investigated by Altena [30], who did not find any 
correlation between membrane properties and the free energy of mixing. For 
this reason the free energy of mixing is not examined in the present work. 


(iit) Viscosity of the casting solutions 

As is generally found [5-8,10,11,13,15,31], an increase in solution viscosity 
due to an increase of polymer concentration hinders the penetration of the 
nonsolvent during the immersion step, and hence the formation of cavities in 
the membrane. Friedrich et al. (10] investigated the effect of the viscosity of 
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Fig. 14. Excess enthalpy of mixing at 37°C of water-solvent systems as a function of solvent molar 
fraction [20]. (See Fig. 13 for symbols. ) 
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equally concentrated polymer solutions in different solvents on membrane 
structure. They showed that the lower is the solution viscosity, the higher is 
the occurrence of cavities in the membrane. For the current system, by com- 
paring the values of the viscosity for the PVDF solutions, listed in columns 4 
and 5 of Table 4, with SEM micrographs of the membranes shown in Figs. 1- 
11, it is possible to check the validity of the conclusions drawn in Ref. [9]. 
When the viscosity value is determined by the polymer concentration, a cor- 
relation holds between viscosity and membrane structure. In fact, with in- 
creasing polymer concentration a trend towards a disappearance of cavities is 
observed (Figs. 10 and 11). However, if different viscosities follow merely from 
the type of casting solvent, membrane structure (Figs. 1-9) appears to be largely 
independent of the viscosity and the results of Friedrich et al. are not confirmed. 


(iv) Diffusivities 

It is likely that in some cases mass transfer could become the controlling 
mechanism of the formation of asymmetric membranes. Unfortunately it is 
not easy to quantify the importance of mass transfer or to model it, since in 
the configurations and systems under investigation mass transfer is a very 
complex process, mainly because of the simultaneous presence of multiple dif- 
fusion, convective transport, thermal effects, nonlinear diffusivities, and phase 
separation. 

Nonetheless it is plausible that, during the very first moments following the 
immersion of the casting solution in the precipitation bath, solvent and non- 
solvent interdiffuse to some extent before polymer diffusion becomes impor- 
tant and/or the skin is formed. The characteristic time for mutual diffusion, 
t,, of solvent and nonsolvent in a typical thickness of the skin (~0.1 um, or 


TABLE 4 


Density and viscosity of PVDF casting solutions 








Solvent PVDF concentration Pes Nes nes 
(%) (kg/m*) (mPa-sec ) 

DMA 15 1008.7 889.1 937.7 
DMF 15 1014.3 639.5 751.4 
DMSO 15 1156.1 3053.1 1394.5 
HMPA 15 1085.0 9074.3 2550.1 
NMP 10 1074.8 447.6 245.2 
NMP 15 1096.2 2239.8 1231.1 
NMP 20 1119.4 8943.4 4918.6 
TEP 15 1129.4 3021.2 1802.4 
TMP 10 1247.6 886.9 403.3 
TMP 15 1265.1 5324.3 2426.1 


TMU 15 1032.0 1687.9 1100.1 
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10~-‘ m) can be evaluated from the order of magnitude of the solvent—nonsol- 
vent diffusivity D~10~° cm?/sec, or 10~° m*/sec [32-34]. It follows that 
t,= (10~-")?/10-°=10~-° sec. The diffusivity for the polymer is at least one 
order of magnitude smaller than 10~° cm?/sec [25]. This implies a diffusion 
time for the polymer in the typical thickness of the skin of much longer than 
t,. The inference is that for many systems, in the very first steps before skin 
formation, the polymer can be considered practically stable, whereas solvent 
and nonsolvent interdiffuse. Similar arguments are the basis of some recently 
developed models [36] where a pseudobinary (solvent and nonsolvent) dif- 
fusion is adopted to describe the early stages of asymmetric membrane for- 
mation, and are also consistent with the observation that the net macroscopic 
polymer movement in such systems can be safely neglected, as the final poly- 
mer distribution suggests [37]. It must be pointed out that such models repro- 
duce with good agreement the experimental data and observations. 

The rate of phase separation is determined by the kinetics and depends on 
the degree of penetration in the demixing gap. In any case, macromolecules 
must rearrange before forming a solid phase. This process cannot be consid- 
ered instantaneous, and the delay between attaining the suitable precipitation 
concentration and the actual growth of a solid skin can be evaluated as being 
some orders of magnitude longer than t, [9]. 

In such a scheme it makes sense to assume that the local concentration at 
formation of the skin, and consequently the skin structure, depends essentially 
on solvent—nonsolvent diffusion. So even if, strictly speaking, diffusion takes 
place in a ternary system, the binary diffusivity of solvent and nonsolvent can 
give some significant insight into the rate of the process. Experimental evalu- 
ation of the binary concentration dependent diffusivities for the different sol- 
vent—nonsolvent systems is difficult and time consuming, but a simple predic- 
tion of an average value should be sufficient for the present purpose and 
qualitative reasoning. As previously illustrated, this prediction was made using 
the Wilke—Chang technique [22]. The computed values at 25°C are reported 
in Table 1 for the different solvents, together with the harmonic mean value 
Dyy=2Ds_wDw_s/(Ds_wt+Dw_s) between Dg w and Dy_p. 

The values of D,,, can be sorted to identify an order of solvents with increas- 
ing mean diffusivity. The resulting order nicely reproduces the same sequence 
as previously detected for decreasing water flux, water content, thickness, or 
cavity volume per unit membrane area. 

A significant dependence of the water flux on the mean diffusivity seems to 
exist, as can be appreciated from Fig. 15, where this quantity is plotted vs. D.,. 
In contrast, plots in Fig. 16 of water flux vs. solubility parameter of solvent, 
[7], mean cloud point, or enthalpy of mixing do not show any firm and precise 
correlation, or look to be largely random. 

Some explanations can be given of the dependence of the resulting structure 
on the solvent—nonsolvent diffusivity. This influence should not be surprising, 














le 
30} 
es 
oe A 
£ 
a 20+ 
x 
Ss 
LL 
A 
10} ° 
O 
_ 1 4 8 EE e 2 
8 10 12 


Diffusivity « 10° (cm?2/s) 


Fig. 15. Water flux vs. mean mutual diffusivity of the system solvent-nonsolvent. (See Fig. 13 for 
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Fig. 16. Water flux vs. solubility parameter of solvent, limiting viscosity number, mean cloud point 
and maximum enthalpy of mixing. (See Fig. 13 for symbols. ) 


as in many models a parameter related to the solvent—nonsolvent diffusivity 
appears. In particular, some models, such as those in Refs. [9] and [36], pre- 
dict that, when the value of the solvent—nonsolvent diffusivity increases, the 
concentration path in the ternary diagram during membrane formation should 
lead to entry into the demixing gap at higher polymer concentration. On the 
other hand, some experimental observations suggest that the occurrence of 
phase separation at high polymer concentration determines the formation of 
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denser and less porous structures than those obtained at lower polymer con- 
centration. By these arguments, the present postulated effects of the solvent- 
nonsolvent diffusivity on the final membrane structure are justified. Nonethe- 
less, as many existing models do not focus particularly on simultaneous mass 
transfer both in the “nonsolvent” bath and in the “casting” solution in the 
very first moments of membrane formation, it might be advisable to check 
further on the influence of the solvent-nonsolvent diffusivity parameter 
through a model that takes special account of these aspects. It is the intention 
of the authors to accomplish this task in future work. 


Conclusions 


The properties and performances of a membrane produced by the phase sep- 
aration technique depend very much on the type of solvent utilized for the 
casting solution. The consequence is that the selection of the proper solvent is 
essential in the design and engineering of a membrane suitable for a certain 
given ultrafiltration process. However, difficulties and uncertainties arise in 
the prediction of membrane characteristics on the basis of the properties of 
the system used to form the membrane. The prediction implies the capability 
of quantifying the relative importance of the mechanisms participating in the 
formation of the membrane, which are mainly phase separation and diffusion 
in the ternary system of polymer-solvent—nonsolvent. Since these two mech- 
anisms are affected to different extents by thermodynamic and transport prop- 
erties, the influence of some of these properties on the membrane structure 
was investigated. To this end the solubility parameters of both polymer-sol- 
vent and solvent—nonsolvent, polymer-solvent compatibility, cloud points, 
solvent—nonsolvent excess enthalpy of mixing, casting solution viscosity, and 
solvent—nonsolvent diffusivity were measured or evaluated. Results indicate 
that, contrary to some other systems investigated in the literature, thermo- 
dynamic properties, examined one by one, do not influence the performance of 
the final PVDF membrane, whereas a good correlation was found between sol- 
vent—nonsolvent mutual diffusivity and the membrane structure. This seems 
to indicate that, in the present system, diffusion is the controlling mechanism 
of membrane formation, in accordance with the conclusions and hypotheses of 
other researchers (see, e.g. Ref. [8] ). 

It must finally be noted that, for PVDF membranes obtained from different 
solvents, many membrane characteristics, such as water flux, water content, 
membrane thickness and void volume per unit area, are correlated one to an- 
other as a consequence of the processes involved in the formation of the mem- 
brane. So any of the previous characteristics, besides quantifying the measured 
property, may furnish information on the resulting membrane structure. 
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Abstract 


The integration of continuous fermentation and membrane separation in the same unit to en- 
hance the performance of a bioreactor has been studied. Pervaporation with ethanol selective 
silicone rubber hollow fibers was considered. This is called the ““Continuous Membrane Fermentor 
Separator (CMFS)”. In order to observe the membrane effects, the continuous fermentation ex- 
periments were conducted through the cyclic operations of dead membrane phase (conventional 
continuous fermentation) and live membrane phase (CMFS) under the same geometrical and 
physicochemical conditions. Relative to conventional continuous fermentation, the performance 
of CMFS results in high yeast cell densities, reduction of ethanol inhibition, longer residence time 
of substrate, more glucose consumption, and recovery of clean and concentrated ethanol; these 
are called ‘“‘membrane effects.” A 10-20% increase in ethanol productivity was achieved and per- 
vaporated ethanol concentration was over 7 times higher than that in the fermentor. A mathe- 
matical model was developed and used to determine the effects of design and operation parameters 
of the CMFS. Computer simulation results predict significant improvement in ethanol productiv- 
ity as well as in glucose utilization when the dimensionless membrane volume was increased. 


Keywords: biotechnology; membrane reactor; pervaporation; ethanol fermentation; silicone rub- 
ber membrane; hollow fiber membrane 





Introduction 


Membrane technology may be divided into three categories: membrane 
preparation, separation of various mixtures, and membrane reactor develop- 
ment. The membrane preparation is an essential element of membrane tech- 
nology. In a separation process, the role of the membrane is to act as a selective 
barrier. Once we have the membrane and know its separation performance 
characteristics [1], the next step would be to apply the membrane to a chem- 
ical reaction. This is called the membrane reactor. Membrane application to a 
biochemical reaction is one of the most attractive research areas since the 
membrane reactor using synthetic semipermeable membranes of the appro- 
priate chemical nature and physical configuration can provide numerous ad- 
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vantages to the development of a bioreactor; the biocatalysts (whether micro- 
bial cells or enzymes ) are physically trapped on one side of a suitable membrane, 
the resulting product is recovered by selective permeation through the mem- 
brane, heat sensitive materials in the biosystem can be treated without dena- 
turing, etc. 

The limitation and unstable supply of petroleum resources will direct inter- 
est to the production of fermenting ethanol from renewable agricultural prod- 
ucts or cellulosic wastes in the near future. The fermentation ethanol may be 
used as a liquid fuel or as a feedstock for the production of organic chemicals 
currently manufactured from petroleum [2]. Development in enzymatic or 
acid hydrolysis of cellulosic materials to glucose may supply fermentable car- 
bohydrates at a lower cost [3,4]. 

Ethanol fermentation is a product-inhibited reaction with a maximum tol- 
erable ethanol concentration from 6 to 10%, depending on the microorganisms. 
Traditionally, fermentation ethanol has been produced in batches. Batch fer- 
mentations are inherently inefficient because they remain “down” for a long 
period of time between batches and are quite unproductive during the startup 
period [5]. Moreover, the ethanol productivity is very low (2-2.5 g ethanol/I- 
hr) due to low concentration of yeast cells, end product inhibition, and sub- 
strate depletion. Conducting fermentation in a continuous mode rather than 
in a batch process is advantageous because of higher productivity, better pro- 
cess control, improved yield, etc. However, conventional continuous fermen- 
tors have one major disadvantage. The dilution rate cannot exceed the growth 
rate of the yeast cells, or there will be cell “washout.” Cell density and dilution 
rate in the conventional continuous fermentations were low, thus limiting pro- 
ductivity. Ethanol inhibition is also an important limiting factor because a 
high tolerance of ethanol permits concentrated glucose substrates to be fer- 
mented with higher utilization of glucose. Therefore, improvement in the pro- 
ductivity of end product-inhibited fermentation requires the removal of the 
product as it is formed. And the beer which is produced with an ethanol content 
of 6-10% is conventionally distilled to at least 95% ethanol in order to be suit- 
able as a fuel. In commercial ethanol production, high concentration of clean 
ethanol is desired for economy in the downstream processing. 

There have been several attempts to achieve a simultaneous separation of 
ethanol as it is formed, along with fermentation. Various product removal 
methods have been coupled with fermentation, including vacuum evaporation 
[6,7], flash vaporization [8], solvent extraction [9,10], and membrane per- 
meation. Systems based upon membrane permeation have been suggested in a 
variety of processes, including dialysis [11], ultrafiltration or microfiltration 
[12-14], and pervaporation [15]. Pervaporation has been considered for the 
ethanol-to-water separation with various membranes [16]. The performance 
of silicone rubber has been reported to have a high ethanol-to-water selectivity, 
while most of the other membrane materials exhibit high water-to-ethanol 
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selectivity. Hoover and Hwang [17] have studied the continuous membrane 
column with pervaporation of alcohols and described a mathematical model 
for capillary membrane pervaporation. Permeability data and separation fac- 
tors as a function of alcohol concentration were also obtained and fitted to a 
polynomial equation. The separation factor of ethanol-to-water was 7 to 8 when 
the concentration of ethanol was below 5%. 

When the conventional continuous fermentation is coupled with pervapor- 
ation by a capillary permeator, it is called the CMFS. The objectives of this 
study were to investigate the effects of membrane pervaporation on the con- 
tinuous fermentation process and to demonstrate that the CMFS fermentation 
system would be superior to the conventional continuous fermentation in terms 
of ethanol productivity, utilization of substrate, and downstream processing 
by comparing both processes under the same geometrical and physicochemical 
environment. In addition, a mathematical analysis has been developed which 
adequately describes the CMFS and predicts the performance of the CMFS 
when the system parameters change. 


Theory 


1. CMFS concepts 

The CMFS, illustrated in Fig. 1, has a shell and tube configuration; in other 
words, a bunch of hollow fiber membranes are retained in a cylindrical shell 
container. 

The proposed CMFS consists of two basic segments: (1) a well-mixed fer- 
menting region in the shell side termed as the “fermentor” and (2) a product 
pervaporating region through the lumen of the hollow fibers, termed as the 
“separator.” Fresh substrate is fed to the fermentor. In the fermentor, glucose 
is converted to ethanol. Pervaporation takes place across the interposed mem- 
brane between the fermentor and the separator. The pervaporated permeate 
contains highly concentrated ethanol solution. Ethanol inhibition may be di- 
minished, which can cause faster rates of ethanol production and yeast cell 
growth. When the appropriate amount of residue is taken from the fermentor 
to keep the continuous fermentation operation, there is a mismatch in flow 
rates between fresh substrate and residue; that is, the amount of residue stream 
becomes less than that in conventional continuous fermentation. Comparing 
the same case in conventional continuous fermentation, yeast cells in the fer- 
mentor will be accumulated and the residence time of the substrate in the 
fermenting region will get longer. The degree of flow mismatch depends on the 
design and operation of the CMFS, such as dilution rate, membrane module, 
permselectivity of the membrane, etc. As a consequence, more ethanol can be 
produced by better utilization of glucose because of the reduction of ethanol 
inhibition and higher yeast cell densities. Hence, the CMFS is expected to 
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Fig. 1. Concept of the continuous membrane fermentor separator. 


result in higher yeast cell densities and higher ethanol productivities along 
with the simultaneous recovery of highly concentrated clean ethanol. 


2. CMFS modeling 

The CMFS is the integration of the fermentor and the separator. The fer- 
mentor considered is the continuous fermentor from which the ethanol is re- 
moved by membrane pervaporation. The separator is the capillary permeator 
in which fermentation occurs in the shell side. Combining the information of 
the conventional continuous fermentation and that of capillary permeator pro- 
vides the appropriate expression for the CMFS model. 
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Analysis of a single stage chemostat [18], shown in Fig. 2, will serve as a 
basis for the conventional continuous fermentation. A set of equations relating 
concentrations of yeast cells X, product ethanol E, and limiting substrate glu- 
cose S to the primary independent operating variable, the flow rate of the fresh 
medium, can be derived through the use of material balances as follows: 














dx V; 

_ nee) 1 
dE, vu; 
dt = (Bi— Bo) + vx, (2) 
dS, ov; 
dt =ySi- So) ye (3) 


Here, subscripts i and o represent inlet and outlet conditions, respectively. 
At steady state conditions, the left hand sides of eqns. (1)-—(3) are zero. 
Then, 
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Fig. 2. Schematic diagram of conventional continuous fermentor. 
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values of D, X,, S;, S,, and E, can be obtained by the conventional continuous 
fermentation experiment. 

Hoover and Hwang [17] have described the pervaporation process with the 
overall phenomenological equation, and the capillary permeator has been mod- 
eled in terms of the liquid pervaporation process. Referring to Fig. 3, which 
shows all the different elements at external connections to the CMFS, the 
model will treat the fermentor as the well-mixed region of fermentation broth 
and the separator as the plug flow region of sweep air stream. The flow rates 
and mole fractions of the ethanol E, water W, and sweep air S in the separator 
side are as follows: 


Gr =Gsz +Gpt+Gw (7) 
Ye =Gp/Gr (8) 
Yw =Gw/Gr (9) 


The differential forms of the pervaporation equation for ethanol and water 
were written for a cylindrical geometry for n hollow fibers with radial dimen- 
sion of r, and r; over a differential length dZ [17]. 
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dGw Qw Itn (Pw Xw —Ps Yw/yw) (11) 
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Since membrane permeation is a slow process, the pervaporation process 
can be considered as a pseudo steady state process during the short time period, 
which means the amount of pervaporation remains constant. The amounts of 
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Fig. 3. Continuous membrane fermentor separator model. 
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ethanol and water pervaporated in a short time interval were obtained by in- 
tegrating eqns. (10) and (11) along the length of the separator /. Then, 





dZ (12) 


Ga, <fFen =| Qe 27 n (Pi Xg — Ps Yr/Ye) 
Eo es : In (r,/T;) 





Gis. Gh, dZ (13) 


l 
-| Qw 27 n (Pw Xw —Ps Yw/w) 
od In (r,/r;) 

If the model of the separator is coupled with that of the conventional con- 
tinuous fermentor, the component mass balances of ethanol E, glucose S, and 
yeast cell X in the fermentor side can be modified from eqns. (1)—(3) taking 
into account eqns. (12) and (13) as follows: 














dE 1 

=——_— (Geo — Gg; +v,E,) + VX, (14) 
dt v 
dS, 1 

° =— (v,S; —v,S,)— X. 
dt 5 Wis Uo = Si oO (15) 
dX Uv 

. ed b 4 6 
— (x « , (16) 


Initially inlet and outlet conditions from the steady state of the conventional 
continuous fermentation and inlet conditions of the separator are given. Then 
solving eqns. (14)-—(16) after first computing eqns. (12) and (13) gives the 
dynamic behavior of the CMFS along the given time intervals. These compu- 
tations are repeated until the concentrations of yeast cells, ethanol and glucose 
reach the steady state. At each time increment, the values of component 
permeabilities and any other variables must be reevaluated as necessary. Va- 
por-liquid equilibrium behavior was recovered from the Antoine vapor pres- 
sure equation and activity coefficients from the Wilson equation to the extent 
that these equations approximate the actual data. The program was solved on 


a Wilbur computer using the differential equation solver DGEAR developed by 
IMSL. 


Experiments 


After the membrane module was constructed, a series of permeation tests 
were conducted to verify the permeability values of ethanol and water through 
silicone rubber membrane. They were very close to the reported values by the 
previous investigators [15,17]. 
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1. Yeast strain and culture 

The yeast organism used in the ethanol fermentation studies was Saccha- 
romyces cerevisiae NRRL-Y-2034 obtained from the U.S. Department of Ag- 
riculture’s Northern Regional Research Laboratory (Peoria, IL). For the prep- 
aration of YM broth to cultivate yeasts, the medium was composed of 3 g yeast 
extract, 3 g malt extract, 5 g pepton, 10 g glucose, 20 g agar and 1 | of distilled 
water without adjusting the pH. Every month new slants of the yeast were 
cultured to insure high viability and to lower the changes of strain degradation. 
The composition of the substrate medium was adapted from Cysewski and 
Wilke [7] and is given in Table 1. Ergosterol was supplemented to eliminate 
the oxygen need to retain the yeast cell viability during the anaerobic condi- 
tions [19]. Citric acid and sodium were added as a buffer to maintain the pH 
of 4.0—4.5 for minimizing the chances of contamination. The medium compo- 
nents were mixed and sterilized together at 121°C for 20 min in an autoclave. 


2. Design and operation of the CMFS fermentation system 

Figure 4 illustrates the design of a laboratory-scale CMFS that was fabri- 
cated in such a manner that it could be readily assembled and disassembled 
between experiments in order to clean and sterilize the various components 
and to change its configuration. The CMFS has a shell and tube configuration 
which consists of an outer shell made of 1 in. glass pipe, 100 cm in height. The 
glass shell is joined at the top to a stainless steel Swagelok cross. One side of 
this cross acted as the access to the fermentor. The other side of the cross 
allowed for venting gases from the top of the CMFS. The bottom of the glass 
shell was joined to a stainless steel Swagelok T connector. The side port of the 
T was the exit from the fermentor. The hollow fiber bundle was made of 103 
fibers with an active length of 98 cm and total surface area, based on log mean 
radius, of 1430 cm’. Both end portions of the hollow fibers were potted within 
the 3/4 in. O.D., 1/25 in. thick wall, and the 10 cm long glass tubing end pieces, 


TABLE 1 


Base fermentation media 








Component* g/l 

Glucose 100.0 

Yeast extract 8.50 

NH,Cl 1.32 
MgSO,-7H.O 0.11 

CaCl, 0.06 
Anti-foam 0.2 ml 
Ergosterol 10.0 mg 

Tap water Make up to 11 





*All salts and glucose are reagent grade. 
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Fig. 4. Details of the continuous membrane fermentor separator unit. 


respectively. The membrane used in this study was Silastic Medical Grade 
silicone rubber capillary tubing obtained from Dow Corning Corporation, Mid- 
land, MI. Table 2 lists some information about the tubing dimensions. Dow 
Corning 730 RT'V fluorosilicone sealant was used as a potting material for the 
membranes, because of its strong fuel and solvent resistance. Two bore-through 
stainless steel Swagelok reducing unions (from 1 in. to3/4 in.) for the entrance 
of a 3/4 in. glass tube were connected to the top of the cross and to the bottom 
of the T respectively. The top and bottom of the CMFS were ended with stain- 
less steel Swagelok reducers. The hollow fiber bundle connected with two 3/4 
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TABLE 2 





Silastic tubing dimensions (cm) 





Inside radius 
Outside radius 
Wall thickness 
Log mean radius 


0.015 
0.032 
0.017 
0.0224 
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Fig. 5. Continuous membrane fermentor separator fermentation system. 


1. Fermentor 

2. Feed reservoir 

3. Nitrogen cylinder 
4. Nitrogen purifier 

5. Intermediate tank 
6. Residue receiver 


8. Separator 
9. Condenser 
10. Product receiver 
11. Refrigerated bath 
12. Manometer 
13. Digital thermometer 


7. Gas discharging bottle 14. Rotameter 


P-1. Residue pump 
P-2. Recirculating pump 


P-3. Feed pump 
P-4. Air circulating pump 
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in. glass tubes can be easily charged inside the 1 in. glass shell. All the ferrules 
were made of Teflon. 

The laboratory-scale CMFS fermentation system was constructed to dem- 
onstrate the feasibility of the CMFS concept and to investigate the effects of 
operating conditions on the performance of the CMFS. A schematic diagram 
of the apparatus used in this study is shown in Fig. 5. The continuous fermen- 
tation experiments were conducted through the cyclic operations of dead and 
live membrane phases at various dilution rates by changing the flow rates of 
fresh feed stream. 

In order to start the dead membrane phase, the valves connecting into and 
out of the separator were closed. Setting the fresh feed flow rate, the residue 
stream was controlled to maintain the liquid level in the fermentor. The system 
was allowed to come to a steady state. Samplings of residue were made every 4 
hr to confirm whether or not steady state had been achieved. When the con- 
centrations of yeast cells, glucose and ethanol remained constant over a 12-hr 
period, steady state operation was assumed. At the end of a certain steady state 
of the dead membrane phase, the live membrane phase was started by opening 
the valves into and out of the separator. Then the air circulating pump in the 
separation part was started. As the fermentor level dropped due to the perva- 
poration of ethanol, the residue stream was adjusted to maintain the fermentor 
level for keeping the working volume constant. If the amount of permeated 
ethanol and residual ethanol was matched to the amount of ethanol produced 
in the fermentor, a new steady state was established. This was the completion 
of one cyclic operation of dead and live membrane phases. The dilution rates 
of 0.071 to 0.262 based upon the total fermentor working volume were covered 
by changing the feed flow rate from 30 to 110 ml/hr. The fermenting broth was 
circulated with a reflux ratio >200. Cyclic operation was repeated for every 
dilution rate. 


3. Analysis 

The yeast cell concentration was determined from the optical density mea- 
surement, using a Bausch & Rohm Spectronic 20 Spectrophotometer at a 
wavelength of 520 nm. The calibration curve was prepared by measuring ab- 
sorbance of samples of varying yeast cell concentrations. The yeast cell con- 
centration of the unknown samples was then obtained by measuring the optical 
density and reading the yeast cell concentration directly from the calibration 
curve. Yeast cell viability was determined periodically, using the methylene 
blue staining technique [20]. Concentrations of glucose and ethanol were de- 
termined by liquid chromatography. The instrument was a Hewlett-Packard 
model 1812A liquid chromatograph equipped with a refractive index detector 
and fitted with 700-CH beer analysis carbohydrate column (300 mm x 6.5 mm), 
purchased from Alltech Inc. The column consisted of a sulfonated polystyrene 
resin as the calcium form. HPLC grade water was used as the mobile phase. 
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The column temperature was 90°C, the column pressure was 40 bar, and the 
eluent flow rate was 0.5 ml/min. Calibration curves for glucose and ethanol 
were prepared, using standard solutions. 


Results and discussion 


1. Fermentation kinetics 

The ethanol fermentation kinetics is discussed based upon the results of 
conventional continuous fermentation (dead membrane phase ). Specific rates 
of yeast cell growth yu, and ethanol production v were obtained using eqns. (4) 
and (5). They are plotted against ethanol concentration in Figs. 6 and 7. Two 
scattered points in Figs. 6 and 7 were obtained at very low dilution rates at 
which the concentration of glucose was less than 1 g/|. Specific rates of yeast 
cell growth and ethanol production show strong dependence upon the concen- 
tration of glucose in these situations. Experimental data were fitted to the 
generalized nonlinear fermentation kinetic model described by Luong [21]. 
This model describes the kinetic pattern of ethanol inhibition in terms of the 
following eight parameters. These values were computed by the nonlinear 
regression procedure using ZXSSQ developed by IMSL as follows. 


[, =0.284 Ey =69 a =1.48 Ks =0.18 
Y, =1.23 Ey =70 £B =2.69 Kg, =0.26 


In general, as the ethanol concentration is increased, the specific rates of 
yeast cell growth and ethanol production decrease. At the fermentation tem- 
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Fig. 6. Specific yeast cell growth rate vs. concentration of ethanol obtained from continuous fer- 
mentation experiments at 25°C. 
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Fig. 7. Specific ethanol production rate vs. concentration of ethanol obtained from continuous 
fermentation experiments at 25°C. 


perature of 25°C, ethanol had a noticeable effect on the growth rate at concen- 
trations even above 5 g/l, and the maximum ethanol concentration at which 
the cells did grow was ca. 69 g/l. Ethanol-producing capability did not change 
much below 20 g/] and was considered to reach complete inhibition at 70 g/l 
of ethanol concentration if the curves are extrapolated. Substrate inhibition 
was not considered in this study since it was not encountered, and because the 
glucose concentration was not high and fermentation took place in a well- 
mixed state. However, if much greater concentrations of glucose (greater than 
150 g/l) appeared in the fermentor, the kinetic model may need to treat sub- 
strate inhibition by adding an additional term [22,23]. 


2. Performance of the CMFS 

The CMFS unit was operated not only at the dead membrane phase but also 
at the live membrane phase with pervaporation. The improvement of the CMFS 
fermentation system can be clearly observed since both cases were compared 
under the same geometrical and physical conditions. The goal of this work was 
to find out the membrane effects of the CMFS compared with the conventional 
continuous fermentation. 

Cyclic operations between the dead membrane phase and the live membrane 
phase experienced the stepwise changes of steady state along with the dilution 
rate increase. One steady state shifted very slowly to another steady state. 
Approximately 3 to 4 fermentor volumes of fresh feed supply were needed to 
get to the steady state. During the live membrane phase, flow mismatch be- 
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tween the fresh feed and residue stream was encountered in the range of 1.3- 
1.9 ml/hr, which equals the amount of pervaporation through the membrane. 
Pervaporated ethanol in the live membrane phase is glucose-free, clear, and 
highly concentrated up to 30%, so this stream should be considered as more 
than a product of the conventional fermentor alone. Silicone rubber mem- 
branes have been found to be very durable. There was no tendency for the yeast 
cells to form a layer on the membrane. A low pH of 4.0-4.5 was maintained. 
Yeast viability remained over 95% for the duration of the 19-day experiment. 

Figure 8 shows the results of continuous fermentations with a 10% synthetic 
glucose medium. The glucose, ethanol, and yeast cell mass concentrations are 
plotted against the dilution rate. The solid line shows the result of conven- 
tional continuous fermentation and the broken line represents the CMFS fer- 
mentation system. As would be expected, ethanol and yeast cell concentrations 
decrease with the increasing dilution rates, and glucose concentration in- 
creases accordingly in both cases. However, the CMFS results relative to the 
conventional continuous fermentation show that the yeast cell mass concen- 
tration in the fermentor was increased, highly concentrated ethanol was per- 
vaporated through the membrane, ethanol concentration was kept at almost 
the same level, and more glucose was consumed. Table 3 contains experimental 
data from the continuous fermentation experiments. 

The yeast cell concentration became higher since the cells were accumulated 
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TABLE 3 


Experimental results of continuous fermentation 

















Operating Flow Dilution Concentration Permeate EtOH 
mode* rate rate in residue produced 

feed 

EtOH Glu. Cell EtOH Flow 
conc. rate 

v D E S X P v E 
D 30 0.0714 48.0 0.3 8.32 3.43 
D 41 0.0976 48.0 0.7 8.13 4.69 
D 56 0.1333 45.1 11.2 7.30 6.00 
L 56 0.1333 42.0 0.3 8.80 300 1.87 6.75 
D 65 0.1548 40.2 19.3 6.90 6.19 
L 65 0.1548 40.1 4.1 8.10 294 1.80 7.28 
D 77 0.1833 34.2 30.4 6.02 6.23 
L 77 0.1833 33.3 19.3 7.20 2.59 1.69 7.00 
D 93 0.2214 25.0 49.2 4.90 5.54 
L 93 0.2214 25.0 47.4 5.25 218 1.58 6.26 
D 110 0.2619 12.0 74.2 2.60 3.14 
L 110 0.2619 11.1 73.5 2.70 144 1.30 3.29 





Feed Conditions: E;=0 g/l, S;=100 g/l, X;=0 g/l. 
Fermentor working volume: 420 ml. 
All data are reported at steady state. 
“D: Dead membrane phase. 
L: Live membrane phase. 


in the fermentor due to the flow mismatch; that is, the amount of residue (out- 
let stream) was less than that of fresh feed. It should also be noted that the 
residence time of the substrate became longer as a result of flow mismatch. 
When higher cell density in the fermentor is obtained, the ethanol production 
rate increases. Then, if the ethanol is not removed, accumulation of ethanol 
takes place and its concentration will be increased. Therefore, the fact that the 
ethanol concentration is kept at almost the same level in the CMFS implies 
that the reduction of potential ethanol inhibition is achieved. It should be men- 
tioned that the ethanol permselective property of the silicone rubber mem- 
brane facilitates the reduction of ethanol inhibition (otherwise ethanol accu- 
mulation might occur) and gives highly concentrated, clean ethanol by the 
pervaporation. With an increase of yeast cell densities and the reduction of the 
ethanol inhibition level, along with longer residence time, more consumption 
of glucose can be accomplished. All of these phenomena are called “membrane 
effects,” which are graphically shown as the gap between the solid line and the 
broken line in Fig. 8. 

At the dilution rate of 0.15 hr~', for example, only 82% of the initial glucose 
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was fermented. In other words, it was not possible to ferment a 10% glucose 
feed completely in the conventional continuous culture. In the case of the 
CMFS, glucose consumption was increased to 97% along with higher yeast cell 
density and almost the same ethanol level as compared with the conventional 
continuous fermentation. Permeate ethanol concentration was 290 g/l, which 
was 7 times higher than that in the fermentor. One of the important factors to 
consider is the final product concentration leaving the fermentor and entering 
the recovery system, since the cost of recovery is proportional to the amount 
of liquid to be handled and is inversely proportional to product concentration. 

It can be seen that lowering the dilution rate allows the complete consump- 
tion of glucose where glucose concentration in the fermentor was held near 
zero. The dilution rate for the complete utilization of glucose in the feed was 
shifted from 0.09 in the conventional continuous fermentation to 0.13 in the 
CMFS. The percentage utilization of glucose may be counted in the productiv- 
ity consideration, if necessary. 

The ethanol productivity is shown in Fig. 9 as a function of the dilution rate. 
The ethanol productivity was the sum of the net weights of ethanol in the 
permeate and in the residue divided by the fermentor volume per unit time. 
The solid line is the result of the conventional continuous fermentation and 
the broken line represents that of the CMFS. Although the ethanol concentra- 
tion in the residue stream shown in Fig. 8 decreases with increasing dilution 
rate, ethanol productivity actually increases up to the dilution rate of 0.17 since 
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Fig. 9. EtOH productivity in continuous fermentation. (———) conventional continuous fermen- 
tation; (— — — ) continuous membrane fermentor separator. 
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the productivity is generally the product of ethanol concentration and the di- 
lution rate. After that maximum point, it begins to decrease because, due to 
yeast cell washout, the increase of the dilution rate cannot catch up with the 
decrease of the product concentration. The maximum ethanol productivity of 
the CMFS was 7.3 g ethanol/1-hr, which was a 18% increase over the maximum 
productivity of 6.2 g ethanol/l-hr in conventional continuous fermentation. 
Interestingly, maximum productivity does not necessarily occur at a dilution 
rate corresponding to the maximum conversion of substrate or the maximum 
permeate ethanol concentration. Therefore, optimization compromises may be 
made, taking into account productivity, permeate ethanol concentration, re- 
sidual substrate, etc. 

Potential advantages of the CMFS include increased ethanol productivity 
and the ability to obtain a highly concentrated clean product. The former al- 
lows a more efficient use of fermentor volume and reduces the size of auxiliary 
equipment, such as media sterilizers, pumps, etc. In addition, the removal and 
isolation of the product in the separator could lead to savings in the down- 
stream processing steps since the pervaporated product is a highly concen- 
trated clean ethanol. Furthermore, the hollow fiber module may permit more 
compact and simpler process equipment than any other membrane modules. 


3. Simulation analysis 

The data obtained from the laboratory-scale CMFS fermentation system 
were used to refine and demonstrate the validity of the CMFS model described 
earlier. The points on Fig. 10 are concentrations of ethanol, glucose and yeast 
cells against dilution rates from the CMFS fermentation experiments for which 
the simulation was performed. Computed results of the CMFS model are plot- 
ted with the solid lines. The simulations show that the quantity and composi- 
tion are in good agreement with the experimental results. 

The model of the CMFS could predict the performance of the CMFS fer- 
mentation system when the system parameters change. Volume ratio of sepa- 
rator fermentor (dimensionless membrane volume ) and glucose concentration 
of the fresh feed were chosen for the demonstration, taking into account pos- 
sible membrane effects on productivity improvement and glucose utilization, 
etc. Here, the dimensionless membrane volume was introduced to indicate the 
relative membrane area per fermentor volume. As this quantity becomes larger 
the membrane area grows proportionally for a given fermentor. 

The CMFS performances shown in Fig. 11 were evaluated at four different 
dilution rates by varying the dimensionless membrane volume in the range 0- 
1.0, while the other operating conditions remained the same as those in the 
CMFS experiments. At a certain dilution rate, e.g. following line Number 2, 
ethanol productivity increases as the dimensionless membrane volume in- 
creases due to membrane effects. However, there exists a maximum ethanol 
productivity at which the complete utilization of glucose occurs. In order to 
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Fig. 10. Comparison of the continuous membrane fermentor separator performance between ex- 
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Fig. 11. Influence of membrane capacity on EtOH productivity at various dilution rates. 
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attain this maximum ethanol productivity, a minimum dimensionless mem- 
brane volume is required. This may be a design criterion, since ethanol pro- 
ductivity cannot be improved even though the dimensionless membrane vol- 
ume is increased. In any case, it is indicated that complete utilization of the 
glucose can be achieved by increasing the dimensionless membrane volume. If 
the dimensionless membrane volume is limited at a certain level for some rea- 
son, e.g. 0.1, the dilution rate necessary to reach the highest ethanol produc- 
tivity is 0.18, whereas this value is 0.22 for a dimensionless membrane volume 
of 0.3. This may also be another design criterion. At the limit of the dimen- 
sionless membrane volume tending to zero, we obtain the conventional contin- 
uous fermentation result. 

Figure 12 shows the CMFS performance for 20% glucose feed at the dilution 
rate of 0.15. The dimensionless membrane volume was varied in the range 0- 
1.0 while the other operating conditions remained the same as those in the 
CMFS experiments. A dimensionless membrane volume of 0.6 will be needed 
for the total consumption of substrate. There also exists the maximum ethanol 
productivity, along with the minimum dimensionless membrane volume re- 
quired, which may again be a design criterion. The higher the glucose concen- 
tration in the feed, the greater the improvement of the ethanol productivity as 
the dimensionless membrane volume increases. This is because ethanol inhi- 
bition becomes a major limiting factor as the glucose concentration in the feed 
is increased. 

The demonstration of significant improvement in ethanol productivity, as 
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well as the minimum dimensionless membrane volume for the complete utili- 
zation of glucose will be of major importance. The CMFS concept may allow 
complete fermentation of high substrate concentration even by severely inhib- 
ited microorganisms, assuming that there is a suitable membrane which can 
be fabricated into a well-engineered membrane module. In summary, the mem- 
brane effects strongly depend on the design of the membrane module and the 
CMFS can be optimized to maintain high ethanol productivity and/or good 
glucose utilization, etc., depending upon the desired situation. 


Conclusions 


The effect of ethanol inhibition on yeast cell growth and ethanol production 
by Saccharomyces cerevisiae NRRL Y-2034 was studied by conventional con- 
tinuous fermentation experiments. Both the yeast cell growth rate and ethanol 
production rate decrease as the concentration of ethanol increases. The pa- 
rameters of the generalized nonlinear kinetic relationship model were deter- 
mined, using the nonlinear regression method. At 25°C, maximum specific 
rates of yeast cell growth and ethanol production were 0.4 hr~' and 1.2 hr~’, 
respectively, with complete ethanol inhibition around 70 g/1. 

A laboratory-scale continuous membrane fermentor separator unit was de- 
signed so that it could be assembled and disassembled easily. A silicone rubber 
hollow fiber bundle was used for the pervaporation of ethanol. A medium of 
10% synthetic glucose was fermented in the dilution rate range of 0.071—0.262 
hr~'. Comparing the CMFS fermentation system relative to the conventional 
continuous fermentation, it was observed that ethanol inhibition was reduced, 
yeast cell concentration was increased, residence time of substrate was longer, 
and a highly concentrated clean product was recovered. These are called 
“membrane effects.” 

In the CMFS fermentation system, ethanol productivity is increased by 18% 
over the conventional continuous fermentation. The dilution rate at which the 
complete utilization of glucose occurs is shifted from 0.9 to 0.13. The permeate 
stream was free from microbial cells and any other substrate, and was a highly 
concentrated ethanol solution, with a concentration up to 300 g/l. The exper- 
iments proved that the CMFS can provide stable performance up to 19 days. 

A mathematical model was developed to determine the kinetics of ethanol 
fermentation and to simulate the CMFS performances satisfactorily. Model 
computations indicate that the improvement in ethanol productivity and/or 
complete utilization of highly concentrated sugar substrate could be achieved 
by a greater membrane capacity. 

This study will give a basis for the rationale for the design and scale up of 
the CMFS for ethanol production and provide a framework for the evaluation 
of the CMFS for application other than ethanol production. Attention needs 
to be focused on the CMFS-based fermentation as a means of improving the 
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productivity of product-inhibited fermentations, such as acetone—butanol- 
ethanol by Clostridium acetobutylicum [25], acetic acid [25], etc. 

For the industrial-scale CMFS fermentation system, it is expected that stan- 
dard size modules would be fabricated and these would be connected in parallel 
to obtain the required capacity for a given design. 
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Abstract 


The appearance of stripes near the membrane in cross-flow filtration of dyed solutes in water 
is modeled as a hydrodynamic stability problem of parallel shear flow in a channel of small width. 

Spanwise harmonic perturbations in thickness of the concentration boundary layer give span- 
wise variations in permeate flux which cause cells of streamwise vorticity that are sustained by 
the shear flow. States of neutral stability are determined by linear theory, giving the spacing of 
stripes as function of a modified Reynolds number based on permeate flux, channel width, and its 
ratio to the boundary layer thickness. The present simple model predicts instability to occur at 
considerably higher values of modified Reynolds number than those prevailing in experiments 
showing stripes. 


Keywords: concentration polarization; theory; ultrafiltration; hydronamic stability 





1. Introduction 


In the process of membrane filtration a liquid solution at elevated pressure 
flows through a channel whose walls are permeable to the solvent. As solvent 
permeates through the membrane the solute builds up as a concentration 
boundary layer (polarization ) along the membrane. Mass transfer and result- 
ing permeate flux depend on resistances in boundary layer and membrane. 

The practice of this type of cross-flow filtration, the design and fabrication 
of membranes and the performance of plate-and-frame filtration units is quite 
well known [1]. The transport processes in fluid and membrane and the foul- 
ing of membranes, however, are not well understood [2,3]. 

For many plate-and-frame applications the flow rate and channel dimen- 
sions (typically 16 1 mm or 160.5 mm) are such that an adequate descrip- 
tion is provided by laminar flow between parallel, permeable plates. Simple 
laminar boundary layer theory predicts a permeate flux that decreases with 
axial position, say as the minus one-third power of distance in the Levéque 
approximation which is appropriate for large Peclet numbers. 

Experimental results [4] for hollow fiber ducts of circular cross section do 
confirm a decreasing flux versus position, but property variations, notably the 
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increase of viscosity with concentration, tend to diminish the decrease. How- 
ever, experimental results for plate-and-frame ducts [1] show nearly constant 
flux versus position. Apparently, the axially increasing concentration polari- 
zation is counteracted by some other phenomenon, such as local mixing. 

Further insight has been gained by visualization experiments using large- 
molecule solutes that are color dyes, such as brown (iron) dextran or blue 
dextran. Employing flow through a narrow rectangular duct with membrane 
on one side-wall, the other side-wall being transparent, patterns of non-uni- 
form coloration have been observed elsewhere [1,5,6,7] and in our laboratory. 
The typical pattern is that of steady colored stripes near the membrane, aligned 
with the direction of the flow (Fig. 1). Such structures suggest convection cells 
similar to those of Gértler vortex flows, to those driven by buoyancy in a hor- 
izontal fluid layer heated from below, or to those that may be at play in some 
electroplating processes [8]. 

Observations at no flow, and a pressure difference yielding permeate flux 
through the membrane, show the formation of an irregular marbled pattern, 
apparently on or very close to the membrane. The appearance of stripes, fol- 
lowing the sudden establishment of appropriate favorable conditions, is a slow 
process, lasting on the order of many seconds. Observable stripes vanish either 
if the flow rate is increased to transitional or turbulent flow or if the driving 
pressure difference for permeate flux is removed. The disappearance of stripes, 
following a sudden establishment of conditions appropriate for no stripes, is a 
rapid process, lasting only seconds or less. 

The foregoing observations suggest that stripe formation is a coupled prob- 
lem of mass transfer and shear flow. Also, a regular spanwise periodic array of 
cells of streamwise vorticity could yield a nearly constant permeate flux by 
slowing or halting streamwise boundary layer growth. 

Since the intensity of coloration of a dyed solution is proportional to the 
integral of concentration of solute over depth it is conjectured that observed 
streamwise stripes represent stable, wave-like variations in thickness and/or 
concentration in the boundary layer along the membrane. Consequently, since 
the local permeate flux varies inversely with boundary layer resistance, say, 
the boundary layer thickness for constant membrane resistance, the observed 
stripe pattern implies a wave-like variation in permeate flux. This variation, 
relative to the uniform permeate flux for undisturbed boundary layer thick- 
ness, implies in effect a periodic inflow and outflow between boundary and bulk 
(Fig. 2). Such perturbations represent cells of streamwise vorticity of alter- 
nating sign, distorting the bulk flow by bringing high momentum fluid to wave 
troughs. This provides a mechanism for sustaining initial perturbations and 
suggests the striping phenomenon to be a stability problem. 

Clearly, in the absence of permeate flow there is no basis for sustaining the 
hypothesized secondary flows, so no stripes should appear. Also, stripes should 
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Fig. 1. Photograph of stripe pattern (vertical lines, indicating also the direction of flow), viewed 
from above the membrane. Double lines are regularly spaced slits (1.0 mm apart) for collecting 
permeate flow under membrane. 

Test conditions: 0.2% blue dextran in water, 160.5 mm straight channel with CA-600 transpar- 
ent membrane, 4p=4 bar, /=0.25 mm, U,,,=0.083 m/sec, 157 mm downstream of inlet. 
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Fig. 2. Model for striping phenomenon, showing perturbed concentration boundary layer and 
resulting secondary flows and bulk flow (schematic). 


disappear in the absence of bulk flow since this flow is responsible for the 
formation of a steady concentration boundary layer. 

Since the phenomenon is observed for solutions of blue dextran, which are 
largely Newtonian, the effect is probably not rheological, i.e. one of shear in- 
duced normal stresses driving secondary flows. The viscosity increases mark- 
edly with concentration, typically by one or two orders of magnitude in the 
boundary layer of blue dextran, but this effect in itself does not seem to provide 
a mechanism for driving secondary flows. 

Presently, the significant increase of viscosity and the small thickness of the 
boundary layer are used in the formulation of a simplified hydrodynamic sta- 
bility model in which cells of streamwise vorticity are driven by boundary layer 
mass transfer and sustained by a channel shear flow. The starting point is a 
simple model for the base flow, treating the mass transfer as a boundary con- 
dition for the bulk flow. The model assumes the base solutions of flow and 
mass transfer to be disturbed by spanwise harmonic perturbations whose sta- 
tionary solutions at neutral stability determine the spacing between stripes. 
Results are compared to observations and limitations of the model are dis- 
cussed in the concluding section. 


2. Physical model and base solution 


The starting point for a stability analysis is to establish the duct flow, con- 
centration boundary layer and permeate flow to be subjected to perturbations. 
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Being concerned with flow in rectangular ducts of large aspect ratio, end effects 
are ignored. 

We therefore consider the steady laminar flow between parallel permeable 
plates spaced a distance 2/ apart. Figure 2 shows the quasi-developed velocity 
and concentration distributions in the symmetric half of this duct flow. Due to 
the pressure difference 4p over the membrane of permeability K, a concentra- 
tion boundary layer of thickness 0, is established. The mass fraction of solute 
is C, and C,,, in the bulk and at the membrane, respectively. The permeate flux 
of solute through the membrane is V,,,. 

To first approximation the base flow is that of steady, incompressible 
Poiseuille flow, 





U=P, (1-7) (1) 
_(=dp/dx)PF 
Po= Qu ] n=y/l, (2) 


where (—dp/dx) denotes the imposed constant pressure gradient. The mean 
velocity U,,,= (2/3)Pp, is of the order 0.1-1.0 m/sec, to be compared to a per- 
meate flux V,,, of the order 10~° m/sec or less. The weak transversal flow as- 
sociated with permeation at the wall may therefore be approximated by 


1 ; 


which satisfies the conservation of mass for quasi-developed Poiseuille flow 
with weak wall suction, hence slowly changing mean velocity. 

Essential features of the mass transfer process may be described by consid- 
ering the somewhat artificial, steady, one-dimensional problem shown in Fig. 
2. The uniform velocity V,,, is the permeate flux which for osmotic pressure 
limited permeation at the membrane is given by 


Vin =K, (Ap — Tim). (4) 


Here, the excess pressure Ap is assumed to exceed the osmotic pressure 7,,, 
which depends on the concentration at the membrane through a virial expan- 
sion such as, for example, for dextran T-70 [4], 


Tm, = (0.6353 C,, +12.64 Ce, +49.55 C2) 10°, (5) 


where units are 7,,(Pa) and C,,,(kg/1) ~ [mass fraction]. The solution to the 
convection—diffusion problem in the layer of thickness 6, (Fig. 2), subject to 
eqn. (5) and bulk concentration C, is 


C,, =Cyexp (Pe Vin/ Vo), (6) 
where V,=K,, 4p denotes the pure water flux and 
Pe= Vo do/D 
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Fig. 3. Concentration at membrane C,,, and dimensionless flux V,,/V>, where V, equals the pure 
water flux, versus permeate Peclet number V,0)/D, the ratio of boundary layer resistance to mem- 
brane resistance. 

Data for T-70 dextran, bulk concentration C,,,=0.002, driving pressure difference 4p: (- - -)1 bar; 
(----) 2 bar; (—) 4 bar. 


denotes the ratio of resistances to mass transfer in boundary layer and mem- 
brane, a water flux Peclet number. 

Solutions to eqns. (4)-—(7) for C).=0.002 and three values of excess pres- 
sures, dp=1, 2 and 4 bar, are shown in Fig. 3 in terms of V,,,/ Vp and C,,, versus 
Pe. Examination of the asymptotic behavior of eqns. (4)—(7) shows, as sug- 
gested by Fig. 3, that Pe-—0 leads to C,,+Cy and V_ approaches a constant value 
which in turn approaches the water flux V, for vanishing Cy. Also, Peco 
shows C,,, to have an upper bound while V,,,/ V, >0 as 1/Pe. The latter behavior 
implies 


Vin ¥D/do, (8) 


being a good approximation for Pe>5 according to Fig. 3. Since our experi- 
ments typically show V,,,/ V) < 0.5, eqn. (8) will be used to describe the typical 
behavior. 

To exemplify orders of magnitude of parameters, let Pex 10, D~ 107" m?/ 
sec and V,~10~° m/sec, then 6) + 10~° m from egn. (7). Of course, d) should 
be determined by solving the convection-diffusion problem of a developing 
concentration boundary layer in shear flow, but such an analysis would give 
comparable values of 6, (see, for example, [4] ). For a stability analysis, how- 
ever, it suffices to describe the essential physics, and in the present analysis d, 
and V,,, are retained as the characterizing parameters for mass transfer. 


3. Perturbation equations 


Considering infinitesimal two-dimensional perturbations in y and z, the re- 
sulting problem of flow and mass transfer is described by 
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U+u',V+v',w' ptp’ do +0',V,,+V',.. (9) 


The linearized equations of motion of perturbations become 





—0 10 
ay t > (10) 





du’ dU. . du’ 1 dp’ 07u' 0?u’ 

—- — + V— = —- — 11 
Ot nial dy oy p Ox E | wey 
dv OV... dv’ 1 dp’ 07v' 82v' 

See g at Waray oe ae 12 
an” dy dy p dy . i dy? | — 
dw’ _ dw’ 1 dp’ 07w' 07w’ 

pt A in A i Oh 13 
Ot a dy p dz | dy? ” a | _ 


where U is given by eqn. (1) and V by eqn. (3). Boundary conditions are, at 
the centerline of symmetry, 


y=0: du'/dy=v' =dw'/dy=0 (14) 
and, at the edge of the concentration boundary layer, y=!—4)—0' = l, 
yxl wu’ =dv'/dy=w'=0 (15) 


The no-slip condition u’ (!) =0 of eqn. (15) is a reasonable first-order approx- 
imation for a vanishingly thin boundary layer, being much more viscous than 
the bulk fluid. Also it is consistent with the base flow assumptions. 


Permeation at the wall, however, can not be ignored. Linearizing eqn. (8) 
gives 


yal: v’ /V.,=—0' /d. (16) 


Finally, we consider the condition of a constant axial pressure gradient acting 
on a quasi-developed flow. From the foregoing no-slip condition the channel 
has effectively a laterally wavy wall represented by the concentration boundary 
layer. The effect on the axial flow will be of order 6’ /I, yielding mean values 
that are larger at troughs and smaller at crests. The corresponding effect in 


velocity gradient at the boundary layer edge is evaluated from eqn. (1), linear- 
ized to 


yal: du' /dy=2 Py 0’ /I?. (17) 


The boundary conditions of eqns. (16) and (17) in effect represent the mass 
transfer problem and its coupling to the bulk flow perturbations. 


According to a usual normal mode analysis, perturbations are assumed in 
the form 


(U' ,0’ w' ,p’ 0’ ,V’ ,,) = (U,0,W,p,6,0 ,, )exp (st +ikz), (18) 
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where k=2z/A denotes the wavenumber in spanwise direction, A being the 
wavelength, i the imaginary unit, and s the amplification factor, generally a 
complex quantity; U, 0, W and p are functions of y only and, they are, along 
with 6 and V,,, the amplitudes of the normal modes. An overbar refers to the 
dimensionless notation below. 

Introducing the notation 


C=ky, ¢, =kl, D=d/dl, o=s/(vk*), Re=V,,1/v 
uw’ =u'Ci/(2ReP,), 0 =u'/V a, w' =w'/Vin; d=0'/do 


where Re denotes the usual Reynolds number associated with suction (per- 
meation ) at a wall, eqns. (10)-—(13), after eliminating W, reduce to 


[D°— (fRe/¢,)D—-1—o]u=—, (20) 
{(D*?—1)[D?— (fRe/f,)D—1—o] +f 'Re/{{}5=0. (21) 


Here, f equals V/V,,, of eqn. (3) andf ' its derivative with respect to n=C€/¢,. 
These terms represent the effect of transversal bulk flow associated with per- 
meation. When these terms can be neglected (see below) eqns. (20) and (21) 
reduce to 


(D?—1-—0)u= — (0, (22) 
(D?—1)(D*—1-—0a)0=0, (23) 


which readily accept closed-form analytical solutions. 
From eqns. (14)-(17), boundary conditions are, at the centerline of 
symmetry, 


¢=0: Du=v=D*0=0, (24) 


(19) 


where eqn. (10), in the form Dv+iw=0, has been used, and at the edge of the 
concentration boundary layer, 


Cw: u=Dv=0, (25) 
v=—9, (26) 
(Re’ /Cj)Du=6d, (27) 

where the modified permeate Reynolds number is defined by 

Re’ =Re 1/69 =(V l/V) (Udo). (28) 


4. Solutions to the stability problem 


Assuming the principle of the exchange of stabilities holds for this problem 
we set a=0 to determine the states of neutral stability. Furthermore, noting 
that V< U, and considering V to be of the same order as perturbations, con- 
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vective terms of the base flow associated with permeation should be of second 
order. 

Therefore, for weak permeation, Re/{, «<1 and Re/¢{j <1, we consider first 
the boundary value problem posed by eqns. (22)-—(27) for which analytical 
solutions are, 


v=A[¢ cosh ¢— (1+(¢, tanh ¢,) sinh ¢] (29) 
u=A[g(¢)—g(¢) cosh C/cosh ¢; J, (30) 


where A is a constant of integration and g(¢) denotes the particular solution 
to eqn. (22), 


g(C) =F(246, tanh ¢,)(C? cosh €—€ sinh €) -; C*® sinh ¢. (31) 


Then, inserting eqns. (29) and (30) into the boundary conditions of eqns. (26) 
and (27), the solution to the eigenvalue problem is obtained by requiring that 
the determinant vanishes for the resulting homogeneous equations in the ar- 
bitrary amplitudes A and 0. This gives 


_ O17 { cosh ¢, — (1+¢, tanh ¢, )sinh Ci} 
= g(¢,) tanh ¢,; —g’ (¢;) 


where a prime denotes differentiation with respect to ¢. Considering the limits 
of eqn. (32) for small and large wavenumber gives the asymptotic solutions 


Re’ (32) 





lim Re’ an lim (Re’ /t) =, (33) 


61-0 2 61 — co 


The solution valid for weak permeation, eqn. (32), is shown in Fig. 4, in 
terms of the dimensionless wavenumber, ¢, =ki/ at neutral stability as function 


/ 


Fig. 4. Dimensionless wavenumber ¢, = k/ at neutral stability versus modified permeate Reynolds 
number, Re’ =(V,,,//v)(l/d,). 
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of modified permeate Reynolds number Re’. The nature of the perturbation 
flow associated with stripes is illustrated in Fig. 5, showing the streamlines in 
the y,z-plane for ¢,=1, 5, 10 and 50. Note that the depth to which the secondary 
flow penetrates into the bulk decreases, as does the wavelength, with increas- 
ing values of wave number. 

Next, including the effect of transversal bulk flow associated with permea- 
tion, eqns. (20) and (21), subject to eqns. (24)—(27), have been solved nu- 
merically. Sample results for two values of the usual permeate Reynolds num- 
ber, Re=0.1 and 0.5, are shown in Fig. 6 in terms of dimensionless wave number 
versus modified permeate Reynolds number. The figure also shows the solu- 
tion for weak permeation, eqn. (32). At low wavenumbers the transversal flow 
has a negligible effect while at high wavenumbers it acts stabilizing, as ex- 
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Fig. 5. Streamlines of perturbation flow in y,z-plane (streamfunction normalized to max. value). 
Cases of dimensionless wavenumber ¢, = 1, 5, 10 and 50. 
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re Z ! 


Fig. 6. Dimensionless wavenumber ¢, at neutral stability versus modified permeate Reynolds num- 
ber, Re’. 

Effect of transversal bulk flow owing to permeation. 

(—) no flow; (----) Re=0.1; (- - -) Re=0.5. 


pected for suction. However, for a value of Re<0.01, which is appropriate for 
the present study, the effect is negligible. 


5. Experimental observations 


Because available experimental results are limited it is only possible to es- 
timate the order of magnitude of the parameters appearing in the stability 
analysis. 

Thus, for a range of flow conditions of a mixture of 0.2% blue dextran in 
water, our observations indicate a typical wavelength of stripes of 0.1-0.15 mm. 
For a channel width of 2/=0.5-1.0 mm this implies a dimensionless wavenum- 
ber of about ¢, ~ 10-30. 

Other quantities known with some certainty include the permeate flux, 
Vn 10~° m/sec and the kinematic viscosity of the bulk fluid v y 10~° m/sec. 
Therefore, the usual permeate Reynolds number is quite small, Re~0.003- 
0.005, and terms of order Re/¢, and Re/{,*, which account for the effect of 
transversal bulk flow associated with permeation, should be negligible. Con- 
sequently, the asymptotic solution, eqn. (33) for ¢, large, should apply, yielding 
the theoretical prediction Re’ ~ (4/3), ~ 13-40. 

Experimentally, 6) is not known directly but estimates based on boundary 
layer theory (see also [4]) suggest d)~0.01-0.05 mm. Hence, the modified 
permeate Reynolds number at which stripes are observed is estimated to be 
Re’ x 0.01-0.25, which is much less than the critical values predicted by the 
model. This discrepancy is probably due to the simplicity of the present model, 
particularly with regard to its lack of detailed description of the mass transfer. 
Also, the model is based on linear, infinitesimal theory. 
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In view of the estimated large values of ¢, the penetration of secondary flow 
into the duct flow will be small (Fig. 5). Considering also the expected mag- 
nitude of these velocity perturbations it is not surprising that they have not 
been detected in exploratory velocity measurements using laser Doppler 
velocimetry. 


6. Discussion and conclusions 


A simplified stability model has been formulated for the striping phenome- 
non observed in the quasi-developed mass transfer process of cross-flow ultra- 
filtration in a parallel plate channel. The boundary layer mass transfer process 
has been simplified, and in effect reduced to boundary conditions imposed on 
the bulk flow. The necessary coupling that sustains perturbations in boundary 
layer thickness, hence in permeate flux, is associated with the primary bulk 
flow which is assumed to be subjected to a constant axial pressure gradient. In 
this way, an increased boundary layer thickness implies a reduced bulk velocity 
and a reduced velocity gradient at the boundary layer. 

The present stability problem, as formulated and solved for the case of weak 
permeation, has some mathematical similarities with that of Gértler vortices 
arising in boundary layer flows over concave walls, see, for example Ref. [9]. 
In the latter problem, the centrifugal body force, through its coupling of uw and 
U, provides a volume source, yielding non-trivial solutions with no distinct min- 
imum on the curve of neutral stability. In the former problem, a volume term 
is absent but the mass transfer provides a coupling between u and U at the 
boundary (eqns. 26 and 27), affording a closed form solution, eqn. (32), also 
having no distinct minimum on the stability curve. 

From experimental observations of stripe spacing and permeate flux this 
simplified case of weak permeation should apply to the experimental condi- 
tions. Hence, neutral stability is governed by a single parameter, the modified 
permeate Reynolds number, eqn. (28). Accordingly, instability is promoted by 
low bulk viscosity, high permeate flux (high permeability and high excess pres- 
sure), and, compared to channel width, a thin concentration boundary layer 
(low mass diffusivity, high bulk flowrate, and low viscosity of wall layer). 

While the foregoing parametric trends appear physically reasonable, nu- 
merical values based on observation (see foregoing section ) do not confirm the 
quantitative results of the model. One reason for the discrepancy may be the 
simplicity of the model formulation with regard to the mass transfer process. 
Another is the linear analysis which may be expected to yield too large critical 
values for the governing parameter [9]. Further study is required to resolve 
the discrepancy. 
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List of symbols 

A integration constant 

C. mass fraction of solute in bulk 

C. mass fraction of solute at membrane 
D mass diffusivity (m*/sec ) 

g particular solution, eqn. (31) 

l imaginary unit, \/—1 

k wave number (1/m) 

K, membrane permeability (m-sec™ '-Pa™ ') 
l duct half-width (m) 

p pressure (Pa) 

Py Coefficient, eqn. (2) (m/sec) 

Ap driving pressure difference (Pa) 

Pe water flux Peclet number, eqn. (7) 
Re _ permeate Reynolds number, eqn. (19) 
Re' modified Reynolds number, eqn. (28) 
s amplification factor (1/sec) 

U,u _ velocity, streamwise (m/sec ) 

U,, bulk (mean) velocity (m/sec) 

V,v __ velocity, transveral (m/sec) 

Vo water flux (m/sec) 

V,, permeate flux (m/sec) 

Ww velocity, spanwise (m/sec ) 

x streamwise coordinate (m) 

y transversal coordinate (m) 

2 spanwise coordinate (m) 


Greek letters 


O boundary layer thickness (m) 
do unperturbed thickness (m) 
n dimensionless coordinate, y/l 


A wavelength, 27/k (m) 
7 dynamic viscosity of bulk fluid (kg/sec-m ) 


py kinematic viscosity, 4/p (m°*/sec) 

4 osmotic pressure (Pa) 

Nm  OSmotic pressure at membrane (Pa) 
p mass density (kg/m) 


a amplification factor, eqn. (19) 








1 





dimensionless coordinate, ky 
dimensionless wavenumber, kl 
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Abstract 


The transport numbers of chloride ion as well as unidirectional fluxes under a constant current 
are measured by the radiotracer technique. The studied anion exchange membranes are in contact 
on both sides with hydrochloric acid at the same concentration. The membranes are the commer- 
cial membranes AAV Selemion and ARA Morgane, specially designed for the recovery of acids by 
electrodialysis. In addition, measurements of the electrical conductances of the membrane and of 
the amount of sorbed electrolyte at equilibrium have been carried out. The analysis of the results 
suggests a low dissociation degree of the acid present in the membrane phase. 


Keywords: ion-exchange membrane; electrochemistry; electrodialysis; Selemion membrane; Mor- 
gane membrane; chloride ion transport; radiotracer technique 





Introduction 


Nowadays electrodialysis can be used for the treatment of acid effluents 
from industry (chemical, hydrometallurgy, electroplating). By this process, 
concentrated acid can be obtained with anion exchange membranes showing a 
low permeability towards protons. 

The recovery of acids by electrodialysis can be performed with sulfuric acid 
[1-5], hydrofluoric acid [6], chromic acid [7] and phosphoric acids [8,9]. 
Concerning hydrochloric acid, some studies have been carried out in order to 
improve the electrodialysis efficiency. For this purpose, bipolar membranes 
[10] as well as modified anion exchange membranes [11] have been used. The 
aim of this work is to compare the transport properties of two polymeric anion 
exchange membranes which have been specially manufactured for the recovery 
of acid by electrodialysis: the AAV Selemion membrane produced by Asahi 
Glass and the ARA membrane produced by Morgane. 

Radiotracers are used for measuring the transmembrane unidirectional fluxes 
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of Cl- when the membranes are maintained between two symmetrical HCl 
solutions. By this technique the values of the Cl~ self diffusion fluxes as well 
as the transport number can be evaluated. The study is completed by mea- 
surements of the membrane conductance and of the amount of the HCl sorbed 
in membranes equilibrated with HC] solutions. The analysis of the experimen- 
tal results provides a better understanding of the factors responsible for the 
reduced permeability of the membranes towards protons. 


Experimental 


Two anion exchange membranes expressly manufactured for the recovery 
of acid by electrodialysis are studied: 

e the ARA membrane, which is an experimental membrane produced by Mor- 
gane and has an exchange capacity of 0.6 mequiv-g~' of dry membrane; 

e the AAV Selemion membrane, a commercial membrane produced by Asahi 
glass, which is 110 wm thick and has an exchange capacity of 0.54 mmol-g~! of 
dry material. 

Before utilization, the membranes are soaked in water for more than 48 hr 
and before the experiments, they are maintained in the studied solution for 24 
hr. Measurements are carried out at 25°C; *°Cl-labelled HCl was obtained from 
Amersham Radiochemical Centre. All the solutions are prepared with deion- 
ized water. 


Self-diffusion flux measurements and transport number measurements 
For these measurements, the same cell has been used; the experimental tech- 
nique has been described elsewhere [12]. 


Electrical resistance of the membrane 

The electrical resistance of the membrane equilibrated with the aqueous so- 
lutions is measured at 1000 Hz by means of a Rhone-—Poulenc clip cell. The 
electrodes are platinized-platinum discs. The solution into which the clip cell 
dips is carefully thermostated. The membrane resistance is the difference be- 


tween the measured values of the resistance with the membrane and without 
it. 


Desorption kinetics and determination of the amount of sorbed electrolyte 

The general principle of the method has been described by Glueckauf and 
Watts [13]. In order to determine the amount of sorbed HCl, the membrane 
is immersed in an aqueous solution of HC] labelled with “°C] until the equilib- 
rium state is reached. The membrane is then removed from the solution and 
its surface is wiped with blotting paper, and the membrane is then immersed 
in a given volume of pure water with a vigorous stirring. By measuring the 
variation of radioactivity in the water at regular intervals of time, it is possible 
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to determine the kinetics of release of the sorbed electrolyte in the membrane 
and its total amount. 


Results 


The water content (WC) is defined by the following formula: 


My — Mp 
ne OE 


where my is the weight of the membrane sample equilibrated with HC] 10~' 
N aqueous solution and mp is the weight of the dried membrane in the chloride 
form, the membrane having been maintained under vacuum at 35°C until a 
constant weight was attained. For the two membranes, the water content cor- 
responds to about eight water molecules per anion exchange group. In Table 1 
are collected the main characteristics of the two studied membranes. 


Exchange ionic fluxes and unidirectional fluxes 

The values of the unidirectional fluxes, at equilibrium and under an applied 
current, are obtained from the slopes of the curves of the variation of the num- 
ber of chloride ions appearing in the second aqueous side, the labelled chloride 
ions having been added to the first aqueous side. An example of such a deter- 
mination is shown in Fig. 1. 

During these experiments, it was not possible to note any significant water 
transport. 

The values of the zero current exchange fluxes (J, ) and of the unidirectional 
fluxes J, under a constant electrical current, are given in Table 2. When a 
constant current is applied, the value of the total flux for a given ion trans- 
ported through the membrane is equal to the difference between the two uni- 
directional fluxes J and MI where J corresponds to the unidirectional flux from 
the cathode to the anode and J to the flux in the opposite direction. In Table 


TABLE 1 


Properties of the two studied membranes 








Code Type Thickness Water Exchange 
(4m ) content capacity 
(%) (mequiv-g~') 

AAV Selemion Anion 114 8.6 0.54 

exchange 

membrane 
ARA Morgane Anion 159 8.7 0.60 

exchange 


membrane 

















—— + + +> 
50 t(min) 


Fig. 1. Variation of the number of Cl~ ions crossing the ARA membrane equilibrated with two 2 
M HCI solutions as a function of time. Current density J=30 mA-cm ~?. J, is the self-diffusion 
flux; J is the unidirectional flux from the cathode to the anode; J is the unidirectional flux from 
the anode to the cathode. 


2, I is the density of the applied current. In order to compare the behaviour of 
Cl- and SO{- , some measurements have been carried out with H.SO, solu- 
tions; in Table 3 the results obtained for 4 N solutions are given. 


Membrane conductivity measurements 

The electrical conductivities of membranes equilibrated with different HCl 
solutions have been measured. In Table 4, the values of these conductivities 
are given for various external electrolyte concentrations. 


Determination of the amount of sorbed electrolyte 

Figure 2 shows the variation with time of the number of Cl~ anions which 
are released into a given volume of 4 N HC] solution. The membrane has been 
previously equilibrated with a *°Cl-labelled 2 N solution of HCl. One may con- 
sider that, when the plateau is reached, the total amount of Cl~ initially pres- 
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TABLE 2 


Zero current exchange fluxes (J,) and unidirectional fluxes J and J of Cl~ ions through the 
membranes separating two HC] solutions at the same concentration. J is the unidirectional flux 
from the cathode to the anode and J the flux in the opposite direction 








Cyc) (mol- I=0 1=30 mA-cm~* 
dm~*) 10° Jy . . 
(mol-sec~ '-cm~*) 10° J 10° J 
(mol-sec~'-cm~*) (mol-sec~ '-cm~?) 





AAV membrane 





0.5 2.1 28.0 4.3 
l 1.4 26.8 7.9 
2 2.2 22.5 8.2 
4 3.5 16.4 9.1 
5 1.4 14.6 9.6 
6 1.5 12.8 8.7 
ARA membrane 

0.5 6.1 28.6 3.2 
l 5.0 26.8 2.9 
2 6.4 25.6 17.2 
3 7.3 21.3 28.6 
4 6.7 17.3 26.2 
5 4.6 14.6 22.6 
6 4.8 12.5 51.4 
TABLE 3 


Comparison between Cl~ and SO{~ fluxes obtained with 4 N HC] solutions and 4 N H,SO, 
solutions; current density = 30 mA-cm~* 








Cl- SO; 
AAV membrane 10° J= 16.4 mol-sec™ '-cm~* 10° J=8.5 mol-sec™'-cm~* 
10° J~ 0.9 mol-sec™'-cm~* 10° J~0 mol-sec™'-cm~? 
ARA membrane 10° J= 17.3 mol-sec™'-cm~* 10° J=9.1 mol-sec™'-cm~* 
10° J= 2.6 mol-sec™'-cm~* 10°J~0 mol-sec™'-cm~? 





ent inside the membrane sample has been transferred into the aqueous phase. 
The total number of Cl~ ions, N%,-, present inside the membrane is related 
to the exchange capacity through the equation: 


C1- =X+N, 


where N., is the number of HCI molecules which are sorbed inside the mem- 
brane. Table 5 shows, for the two studied membranes, the values of N %,- and 
N, expressed as mequiv-g~' and mmol-g~' of dry membranes. 
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TABLE 4 


Electrical conductivity of membranes equilibrated with HC] solutions: K x 10° Q-'-cm~' 











[HCl] AAV ARA 
(mol-1~') membrane membrane 
0.5 2.97 7.03 
1 3.50 7.50 
2 4.02 9.19 
3 4.68 11.19 
3.5 4.98 11.53 
4 2.45 11.19 
5 2.76 10.75 
6 3.23 9.75 
nx 104 











. sis , : os - 
50 t(min) 


Fig. 2. Variation with time of the number of Cl~ ions released into 4 N HCI solution from an ARA 
membrane previously equilibrated with 2 NH*°Cl solution. 





TABLE 5 


Amount of total Cl~ ions of sorbed HCl in the membranes 

















| HCl] AAV membrane ARA membrane 
(mol-l~*) 

N* N, N* N, 

(mequiv-g~') (mmol-g~') (mequiv-g~ ') (mmol-g~') 
0.1 1.18 0.64 0.74 0.14 
0.5 1.50 0.96 0.91 0.31 
l 1.94 1.40 1.22 0.62 
2 2.18 1.64 1.33 0.73 
3 2.55 2.01 1.43 0.83 
4 2.67 2.13 1.66 1.06 
5 3.21 2.67 1.95 1.35 
6 4.54 4.00 2.48 1.88 
Discussion 


When aconstant current is applied, the value of the total flux for a given ion 
transferred through the membrane is equal to the difference between the two 
unidirectional fluxes J and J. The transport number of an ion 7 through the 
membrane is then given by the following equation. 


> « 


- |J;—d; |F 


t; I 





(1) 


where / is the density of the applied current and F is the Faraday constant. 
The values of t,;- and tgo3- deduced from eqn. (1) are given in Table 6. 


TABLE 6 


Values of Cl~ and SO{~ transport numbers 








AAV membrane ARA membrane 

[HCl] (mol-1~') 

0.5 0.92 0.94 
l 0.87 0.87 
2 0.74 0.82 
3 0.62 0.66 
4 0.54 0.55 
5 0.46 0.48 
6 0.41 0.41 
|H,SO,] (mol-l~') 

2 0.28 0.59 
4 0.16 0.36 
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In order to analyse the results of transport number measurements, the gel 
model for the membrane phase can be used. The membrane is assimilated to a 
polyelectrolyte solution; considering that the moving electrical species are only 
H,0* and Cl~, the expressions for the transport numbers in the case of com- 
plete dissociation are: 


mM, UH30 + 








t fie = = 2 
_ MUy,0+ + (M+X) Uci- (2) 
and 

(X+mM) dc 
tc- — = (3) 


(X+mM) Uc- + Miiy,0+ 


where m is, at the sorption equilibrium, the molality of the electrolyte which 
is present in the membrane, X is the molality of fixed charges, and a is the 
apparent mobility of ions in the membrane. 

Assuming that the ion mobilities do not depend on the molalities of the 
sorbed HCl, eqn. (3) gives: 








l 1y.0+ 
sani at._2 (4) 
to- Ucj- 1+8 

, m Nyc 

Pp Nsites 


The £ values are given in Table 7. 

One notices the very high values of f, especially with the AAV membrane, 
even at low concentrations. In Fig. 3, 1/t.,;- is plotted as a function of [f/ 
(1+ £)]; it may be observed that the experimental results do not fit eqn. (4) 
in the molality range considered. 

This result must be related to the results of conductivity measurements. 


TABLE 7 


Values of 8 for HCI solutions 








| HCl} AAV membrane ARA membrane 
(mol-l~') 

0.1 1.19 0.23 

0.5 1.78 0.52 

l 2.59 1.03 

2 3.04 1.22 

3 3.72 1.38 

4 3.94 1.77 

5 4.94 2.25 

6 7.41 3.13 

















Fig. 3. Variation of 1/t,,- as a function of B/(1+ £). 


Figure 4 shows the variation of the membrane conductivity as a function of the 
external HCl concentration. A sharp maximum is observed for 1.4 with the 
ARA membrane and a smooth maximum is observed for S=3.8 with the AAV 
membrane. 

The results show that the theoretical model —- the membrane considered as 
a gel phase and complete ion dissociation — does not fit to the studied system. 
One notes the very high amount of sorbed electrolyte in the membrane with 
respect to the amount of water. For instance, in the membranes equilibrated 
with 6 N HCI solutions, the number of sorbed HCl molecules over the number 
of fixed sites (N./X) is 3.13 for the ARA membrane and 7.4 for the AAV mem- 
brane; the number of water molecules over the number of total chloride ions is 
3.46 for the ARA membrane and less than 1.9 for the AAV membrane. 

These data strongly suggest an incomplete ionic dissociation inside the 
membrane phase. A fraction of sorbed HCl probably remains in the neutral 
molecule form. Under these conditions, another expression applies: 











(X+ym) Uc- 
(X+ym)Uc- +YMUy,0+ 


(5) 





tc)- = 


where y is the dissociation coefficient of sorbed HCl. Equation (5) can be rear- 
ranged into the following form: 





(6) 


Equation (6) shows that the values of the transport number of the moving 
ions inside the membrane phase are dependent not only on the ratio of their 
mobilities and the amount of sorbed acid but also on the dissociation coeffi- 
cient of sorbed acid. 

If the ratio of the mobilities of the free ions Wyj4,9+/Uc;- maintains a constant 
value, 1/t-,- will have the same variations as yf/(1+ yf). In Fig. 5 the calcu- 


K.10°(a-.cm”) 
A 


20+ 


“aie 
10f | ag 








Fig. 4. Variation of membrane conductivities as a function of f. 
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lated variations of this expression are plotted as a function of # for four differ- 
ent values of y. From this curve, we can deduce that 1/t,)— is lowered when the 
dissociation coefficient y is low. 

We can conclude that the lower dissociation of sorbed HCl is a factor which 
decreases the proton leakage of the anion exchange membrane. 

As for the variation of the ratio of the mobilities of H,0* and Cl~, no direct 
information can be deduced from the experimental data. We can only observe 
a maximum in the curves of the membrane conductivity (Fig. 4). Considering 
again an incomplete dissociation of the sorbed hydrochloric acid in the mem- 
brane phase, we can only conclude that the mobilities of the free ions have 
decreased because the expression of the membrane conductance is: 


A=F[(X+ ym) tg- +yMty,0+ | (7) 


The decrease of the values of free ion mobilities has two distinct origins: the 
increase of the local viscosity of the hydrophilic pathways in the membrane 
phase, and the increase of the interionic interactions, these two effects being 





— 
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Fig. 5. Variation of yB/(1+y) as a function of £ for several values of y. 


»| 
a 





68 


due to the high value reached by the mobility of hydrochloric acid in the mem- 
brane phase. 


From the data given in Table 2, one may note that, for a current density of 
30 mA-cm ~”, the unidirectional flux J of Cl~ through the ARA membrane 
from the anode to the cathode steadily increases as the amount of sorbed elec- 
trolyte increases. This result means that chloride ions are associated with the 
movement of positively charged species. This fact may be due to the formation 
of an aggregate form such as (H,OC1)* resulting from the solvation of a proton 
by a water molecule and an HC! molecule, ion association inside the membrane 


overcoming the state of a neutral HCl molecule. This result confirms the role 
of ion association in this membrane. 
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Abstract 


The coagulation of ethylene vinyl alcohol (EVAL) copolymer membranes was investigated with 
a solvent/nonsolvent mixture. A two step mechanism model is presented which describes mass 
transfer processes occurring for specified membrane formation systems. The various morphology 
structures can be reasonably explained by this model by changing the solvent content of the co- 
agulation bath. It clearly shows that the membrane structure actually depends on the local com- 
position at phase transition and is controlled mainly by kinetic and structural factors rather than 
thermodynamic equilibrium. This model also points out some aspects neglected of membrane 
formation in published research. This analysis provides better understanding and correct control 
of membrane formation techniques. 


Keywords: membrane preparation and structure; two step mechanism model; phase transition; 
ethylene vinyl alcohol membrane 





Introduction 


Since the development of asymmetric type membranes by Loeb and Souri- 
rajan [1], the field of membrane science and technology has experienced an 
increasing growth. The asymmetric membranes are generally produced by the 
so-called phase inversion technique [2]. In this process, a polymer solution is 
immersed into a coagulation bath containing nonsolvent miscible with the sol- 
vent. A solvent/nonsolvent diffusional interchange between the polymer so- 
lution and the coagulation bath causes the polymer to precipitate, and a ‘gel’ 
containing considerable amounts of solvent and nonsolvent is formed. Finally, 
the gel solidifies to form a membrane film. 

From the technological point of view, the most important problem is to ob- 
tain membranes of optimum morphological structure to meet with the require- 
ments for a wide range of applications, such as gas permeation and blood di- 
alysis. In a previous paper [3], we have described the relationship between the 
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structure of EVAL membranes and experimental variables, including temper- 
ature and solvent effects. Our present discussion will be confined to an isolated 
problem: the addition of solvent to the coagulation bath in an attempt to gain 
a clear understanding of the membrane formation mechanism. 


Theory 


It is obvious that an explanation for the phase inversion technique should 
be based on two factors: (1) the thermodynamic equilibrium properties, in- 
volving the ternary phase diagram for the system polymer/solvent/nonsolvent 
(P/S/NS), and (2) the kinetic properties, involving the mass transfer rates 
of solvent and nonsolvent. 

An isothermal ternary phase diagram is represented in Fig. 1. It is divided 
into three regions according to their thermodynamic state. Two different types 
of phase transition can be distinguished. 

(1) Liquid-liquid phase separation (from I to II): the completely miscible so- 
lution crosses the binodal boundary (the solid line, see Fig. 1), i.e. from a stable 
homogeneous solution region into a two-phase region. 

(2) Solidification (from I to III or from II to III): a hypothetical boundary is 
located in the diagram since the viscosity of polymer solution increases to a 
certain value, the polymer motion will be limited, and then the membrane 
structure is specified. This boundary (the dotted line, see Fig. 1) moves across 
the diagram from left to right while descending slightly since the viscosity in 
concentrated polymer solution is higher in a poor solvent than in a good one 
(not the same as suggested by Mulder et al. [4] ). This is rationalized by pos- 
tulating tighter coiling in good solvents as concentration is increased, where, 
since polymer-solvent interaction > polymer-polymer interaction, intermo- 
lecular contact between polymer chains is forbidden. The coils must as a result 
become smaller as they are crowded together with increasing concentration. 
The result is that each coil contributes progressively less to overall viscosity. 
In poor solvents, on the other hand, intermolecular contacts are not so strongly 


Polymer (EVAL) 





Solvent Nonsolvent 
(DMSO) (water) 





Fig. 1. Ternary phase diagram of EVAL-DMSO-WATER system. I: homogeneous solution region; 
II: two-phase region; III: solidification region. (data are given in Ref. [3] ). 











Fig. 2. Schematic route of composition for dense structure (1), and microporous structure (2) of 
a membrane. 


forbidden and the increasing size of the polymer network owing to chain in- 
terpenetration results in increased viscosity [2]. (The above result is opposite 
in diluted polymer solution.) From the phase diagram it follows that the re- 
sulting structure of the membrane does not only depend on the average com- 
position but also on the route by which this composition was reached (Fig. 2). 
Therefore, the morphology is determined by the local polymer solution in each 
layer of nascent membrane gel, whether the coagulation path follows the line 
1 or 2 (Fig. 2). The line 1 (from I to III) occurs at higher polymer concentra- 
tions. A dense structure is formed in this way. Line 2 (from I to II) occurs at 
lower polymer concentrations. The result is a porous structure. Therefore, it 
is important to know how the rates of transfer of the individual mobile com- 
ponents affect the route of composition change. 

A very important controllable factor is the addition of solvent to the coag- 
ulation bath in order to change its composition. This does not effect phase 
equilibrium but does control the kinetics of mass transfer. Generally, an in- 
crease of solvent content in the coagulation bath reduces the driving force for 
diffusion. The same observations have been made by Strathmann et al. [5], 
Wijmans et al. [6] and Yilmaz and McHugh [7]. But their results are either 
restricted to a short distance from the membrane surface or not consistent with 
the overall membrane properties. The aim of this study is to get a better un- 
derstanding of the membrane formation process by providing a two step mech- 
anism model based on experimental observations. 


Experimental 


Materials 

The ethylene-vinyl alcohol (EVAL) copolymer which has a content of 56 
mol% of vinyl alcohol monomeric units was obtained using previously reported 
techniques [3]. De-ionized and ultrafiltrated water was used in our experi- 
ments. Other reagents are of chemical reagent grade and were used without 
purification. 
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Membrane preparation 

EVAL copolymers were dissolved in DMSO to form a 25% polymer solution 
at 60°C. The resulting solution was kept at 25°C for 24 hr and was subse- 
quently spread on glass plates to a uniform thickness. The glass plates were 
immersed immediately in a 25°C coagulation bath of varying DMSO/water 
composition (see Table 1). After standing overnight, the membrane was re- 
moved from the glass plate and was kept in a pure water bath. 


Membrane characterization 

The morphology of the membranes was examined using a scanning electron 
microscope (Philips Model-515). The water permeability was measured using 
a type 8050 Amicon Filter at 25°C. 

The water content of the membrane was determined from the weight differ- 
ence between an equilibrated membrane in a pure water bath and that of the 
membrane dried to constant weight at 60°C. This content equals approxi- 
mately the percent porosity of the membrane and indicates a diffusional inter- 
change between the polymer solution and the coagulation bath. 


Results and discussion 


Porosity and absorption equilibrium 
Figure 3 shows the liquid content of the membranes as a function of the 
DMSO content of the coagulation bath. It shows that the liquid content after 


TABLE 1 


Details of coagulation baths (D: DMSO, W: water) 





Notation of Gl G2 G3 G4 G5 G6 
membrane 
D/W ratio 0/1 1/4 2/3 1/1 2/1 3/1 
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Fig. 3. The effect of the composition of the coagulation bath on the liquid content of EVAL 
membranes. 
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Fig. 4. Absorption equilibrium of dense EVAL membranes in the water/DMSO bath. 


coagulation in a pure water bath (membrane G-1, see Table 1) is more than 
the original liquid content (75%; see the dashed line), which indicates that the 
DMSO loss by diffusion is less than the water gain. As the DMSO content of 
the coagulation bath increases, the liquid content initially decreases slightly, 
then goes through a minimum, and finally rises. As the DMSO content of the 
coagulation bath is 20% (membrane G-2), the liquid content is about equal to 
the original liquid content, which indicates that the diffusional interchange is 
about the same. 

The minimum point is at a DMSO content of the coagulation bath of 66% 
(membrane G-5). The original liquid content is surpassed as the DMSO con- 
tent of the coagulation bath reaches 75%. In fact, the membrane is swollen by 
the coagulation bath and almost dissolves the polymer. Membranes produced 
in this way have very little mechanical strength. 

Figure 4 shows an absorption equilibrium curve as mentioned below. A dry- 
dense EVAL membrane, which is free of any voids by SEM visualization, is 
immersed in coagulation baths of different DMSO concentrations. Because 
there are not many pores in the membrane, the solvent can only slowly dissolve 
into the membrane. Membranes are weighted after equilibration and the change 
in weight from the dry-dense membrane is drawn in Fig. 4. It shows that the 
absorption equilibrium amount increases monotonously with the DMSO con- 
tent in the coagulation bath. 

The outcome of Fig. 4 is different from Fig. 3. No minimum peak can be seen 
in Fig. 4. It shows that membrane formation mechanism and morphology 
structure do not change monotonously according to the solvation power (i.e. 
the DMSO content) of the coagulation bath. 


The chemical potential variation of the coagulation bath 

The SEM pictures (Figs. 5-10) reveal that the DMSO content exerts an 
influence upon the structure of membranes. This can be explained by the fact 
that due to the presence of DMSO the membrane formation mechanism is 
changed. Therefore, it is very important to know the effect of the DMSO con- 














Fig. 5. SEM picture of the G-1 membrane in a 100% water coagulation bath: (a) top, (b) bottom 
and (c) cross section. 





Fig. 6. SEM picture of the G-2 membrane in a water/DMSO coagulation bath, containing 20% 
DMSO: (a) top, (b) bottom and (c) cross section. 


tent on the chemical potential changes of the coagulation bath since it directly 
influences the driving force for the diffusional interchange between the coag- 
ulation bath and polymer solution [6]. The relation between flux J; and chem- 
ical potential difference 41; can be represented by a simplified phenomenolog- 
ical equation: 


J, =L; Au; (1) 


L; is the permeability coefficient which may be a function of concentration. 
Since the polymer solution system is fixed, we only consider the chemical po- 
tential variation of the coagulation bath (we disregard the influence of L;). For 














Fig. 7. SEM picture of the G-3 membrane in a water/DMSO coagulation bath, containing 40% 
DMSO: (a) top, (b) bottom and (c) cross section. 





Fig. 8. SEM picture of the G-4 membrane in a water/DMSO coagulation bath, containing 50% 
DMSO: (a) top, (b) bottom and (c) cross section. 


the chemical potential of DMSO and water as a function of DMSO content in 
the coagulation bath, the expression used is: 


ui; =U? +RT In y; x; (2) 


where 1? is the chemical potential of component i in a pure (standard) state 
at the same temperature and pressure as that of the mixture. The index  in- 
dentifies DMSO or water. The activity coefficients, y;, used to calculate the 
chemical potential, are taken from Lam et al. [8] and listed in Table 2. Using 
these quantities, (u;— “? ) was calculated from eqn. (2) and is shown in Table 
2. Examining the results (see Table 2 or Fig. 11), one finds that the chemical 
potential of DMSO increases very rapidly and that of water decreases slightly 














Fig. 9. SEM picture of the G-5 membrane in a water/DMSO coagulation bath, containing 66% 
DMSO: (a) top, (b) bottom and (c) cross section. 





Fig. 10. SEM picture of the G-6 membrane in a water/DMSO coagulation bath, containing 75% 
DMSO: (a) top, (b) bottom and (c) cross section. 


upon increasing the DMSO concentration. It is clear that the outdiffusion rate 
of DMSO is reduced very rapidly and the indiffusion rate of water is reduced 
slightly by adding solvent to the coagulation bath (Fig. 12). Therefore, the 
DMSO outflow decreases appreciably and the water inflow decreases slightly. 
In other words, the DMSO outflow decreases more rapidly than the decreasing 
water inflow. This result should lead to a more porous polymer structure in the 
resulting membrane with increasing DMSO content in the coagulation bath. 
However, it provides a poor prediction for the membrane porosity (see Fig. 3) 
while presenting a good fit for the determination of skin structure from SEM 
pictures (see Figs. 5-10). So we investigate the overall membrane formation 














TABLE 2 


Activity coefficients and chemical potential of DMSO-water system at 25°C 


























D/W ratio Yo Yw Lp — Mp lw — Hw 
(cal/mole ) (cal/mole ) 

0/1 0.05 1.00 —o 0.0 
1/4 0.10 0.99 — 3029 — 43.8 
2/3 0.16 0.93 — 2249 — 132.3 
1/1 0.24 0.86 — 1810 — 221.4 
2/1 0.49 0.66 — 1061 — 492.0 
3/1 0.65 0.56 — 755 — 676.2 
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Fig. 11. The chemical potential of DMSO and water relative to their pure state in the coagulation 
bath. 


decrease slightly 





increase rapidly Lo" # 


interface 
coagulation bath polymer solution 


Fig. 12. The variation of the chemical potential of DMSO and water with increasing DMSO con- 
tent in the coagulation bath. 


mechanism in the discussion below. This analysis is an important step towards 
a better understanding of the diffusion process during membrane formation. 


The mechanism of diffusion-controlled formation 
Immediately the casting solution becomes turbid when DMSO is absent in 
the coagulation bath, caused by a rapid diffusion of DMSO from the polymer 
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solution to the coagulation bath. At the same time, relatively little water pen- 
etrates since it is a nonsolvent for the EVAL copolymer. This means that the 
local polymer concentration increases and hence the viscosity also increases 
rapidly at the surface. When the viscosity increases to a certain extent, the 
homogeneous polymer solution undergoes phase transition into a gel and the 
system is solidified (line 1 in fig. 2). A high polymer concentration gel is formed 
in this way, i.e. a dense skin is formed. Figure 5 shows the near absence of pores 
in the toplayer of membrane G-1. The skin formed at the surface of the initial 
membrane is believed to control the rate of mass transfer by having a perme- 
ability considerably lower than the inner portions of fresh solution. Thus, the 
structure below the skin is determined by the diffusion rate of solvent and 
nonsolvent through the skin layer. 

The DMSO outdiffusion rate is thus hindered due to the presence of the 
skin. Water, which has smaller molecular size, penetrates the dense skin layer 
by a molecular diffusion mechanism causing phase transition at a lower poly- 
mer concentration (line 2’ in Fig. 2) resulting in the so-called finger structure. 
So the membrane G-1 formation mechanism has two steps and the structure 
is determined by the local composition in each step. In the first step, the DMSO 
outdiffusion rate is faster, and in the second step, water indiffusion rate is 
faster. The skin is very thin, so the liquid content of the membrane G-1 is larger 
than that of the original solution. Its porosity is about 77%. 

With increasing DMSO content in the coagulation bath, the initial tendency 
of DMSO diffusion from the film slows down since the chemical potential of 
DMSO increases very steeply in the coagulation bath. Therefore, the viscosity 
increases less rapidly and it takes more time for the surface solution to undergo 
phase transition and more water will diffuse into the surface solution, thus 
affecting the skin composition. It follows from line 1’ in Fig. 2 that the result- 
ing skin will be less compact. The Figs. 6 and 7 show that some pores occur in 
the top skin layer. 

Below the less dense skin layer, the DMSO outdiffusion rate increases, and 
the finger cavities gradually become smaller as shown in the SEM picture of 
the cross section. Upon increasing the DMSO content to 66% in the coagula- 
tion bath, the finger structure disappears and the sponge structure emerges. 
Membrane G-5 from this coagulation bath does not have a skin layer. Since 
the tendency of DMSO outdiffusion becomes smaller, the phase transition of 
the surface solution is slow rather than rapid and water has enough time to 
diffuse into the polymer solution to cause liquid-liquid phase separation (line 
2 in Fig. 2). The average pore size of the toplayer is about 0.4 wm. Because this 
is larger than the size of diffusing component, the microporous toplayer does 
not influence the diffusion rate, and the solvent and nonsolvent can diffuse 
freely through the toplayer. Therefore, the diffusion rate of solvent and non- 
solvent below the toplayer is the same as that of solvent and nonsolvent in the 
initial period and the polymer solution below the toplayer will go through the 
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same phase separation as the top layer. Therefore, membrane G-5 is regarded 
to have a one step formation mechanism and a homogeneous structure with 
pores uniformly distributed through the entire membrane (see Fig. 9). Hence, 
the DMSO diffusion rate from the interior polymer solution is higher than 
when a skin is present. Finally, the porosity of membrane G-5 becomes smaller 
and is ca. 66%. Upon increasing the DMSO content in the coagulation bath to 
75%, DMSO presents good solvent power. So the porosity of membrane G-6 
increases, no skin and no finger structure appear, but the membrane does not 
have enough mechanical strength for any application. 

It is obvious from the above discussion that with increasing DMSO content 
in the coagulation bath the toplayer structure follows the line 1 to 1’ or 2, and 
the sublayer structure follows the line 2’ to 2 in Fig. 2. This phenomenon can 
clearly be ascribed to the variation of the diffusion rate of solvent and nonsol- 
vent during the membrane formation process. So the membrane formation 
process clearly has two different mechanisms and the porosity as well as func- 
tion of the membranes is determined by both toplayer and sublayer. 


Membrane performances 

Figure 13 shows the water flux experiment at low pressure difference (<0.9 
kg/cm”). The water flux increases with increasing DMSO content in the co- 
agulation bath, except for membrane G-5. This can be explained by the fact 
that the skin becomes less dense with increasing DMSO content in the coag- 
ulation bath. But the porosity of membrane G-5 is lower resulting in a lower 
water flux than for membranes G-3 and G-4. 

Figure 14 shows the dynamic creatinine dialysis experiment with membrane 
G-1, G-5 and the commercial membrane “Gambro”. It is obvious from the 
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Fig. 13. Influence of DMSO content of the coagulation bath on the water flux of EVAL membranes 
at different operating pressure. 
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Fig. 14. The creatinine dialysis experiment for G-1, G-5 and Gambro membrane. 


figure that membrane G-5 is better than Gambro. Therefore, membrane G-5 
is a potential film for the application in dialysis. 


Conclusion 


In this paper we have presented the mass transfer processes responsible for 
the membrane formation mechanism. The ultimate determining factor for 
membrane structure is the local composition of each layer in the casting solu- 
tion at the moment that a phase transition is reached. Progressing solidifica- 
tion actually depends on the diffusional interchange of solvent and nonsolvent. 
As shown, the mass transfer rate influences membrane formation, which can 
be considered as a two step mechanism. 

(1) The ratio of nonsolvent inflow and solvent outflow is of the utmost im- 
portance for the toplayer structure, more specifically, the compactness of the 
toplayer actually depends on the local composition at phase transition. If the 
toplayer enters the solid phase directly, a dense skin is formed. If the toplayer 
crosses the binodal boundary to cause liquid-liquid phase separation, the result 
is a porous layer. The structure of the toplayer is thus related to the chemical 
properties of the polymer/solvent/nonsolvent system, and especially the vari- 
ation in chemical potential has a remarkable influence. 

(2) If the toplayer structure is relatively compact, it will increase the barrier 
for diffusion of solvent from and nonsolvent into the sublayer of polymer so- 
lution. Therefore, mass transfer in the dense toplayer (skin) occurs at a much 
slower rate, and it will determine the overall diffusion rate through the film 
below the toplayer. On the other hand, if the toplayer structure is microporous, 
mass transfer rate is only slightly influenced by the toplayer. As a consequence, 
the sublayer structure is the same as the toplayer structure and the two step 
mechanism is reduced to a one step mechanism. Thus, the sublayer structure 
can be predicted from the chemical and physical properties of the toplayer/ 
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solvent/nonsolvent complex, and especially the toplayer structure plays an 
important role. 


Applying the above result, we are able to govern membrane structure by 
suitably regulating the solvent content in the coagulation bath. Asymmetric 
and symmetric, finger type or sponge type membranes can be obtained by the 
phase inversion technique. Finally, our results demonstrate that the two step 


mechanism model leads to a good prediction, which is consistent with existing 
experimental observations. 
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Abstract 


The responses of an amperometric electrochemical cell to species of interest in the field of 
clinical biosensors, including hydrogen peroxide, ascorbate and phenolics were investigated. The 
effects of covering the electrode surface with polycarbonate and cellulose-based membranes on 
the response were determined, in order to assess the suitability of such membranes for use in 
biosensors. 

The addition of an amino or a further hydroxyl group in the ortho or para positions to a phenol 
group, enhanced the electroactivity of the phenolic group. A polycarbonate membrane of 2 “m 
pore size imparted no selectivity, with all phenolic, ascorbate and hydrogen peroxide electrode 
responses being 10-20% of the bare electrode value. 

Cellulose acetate membranes gave smaller responses, but were selective, being relatively more 
permeable to hydrogen peroxide and to phenol than to other phenolics, and excluding ascorbate 
to a large extent. Such membranes thus have a useful selectivity for use in certain biosensors. 
Cellulose acetate butyrate membranes showed a further reduction in signal size, but an unchanged 
pattern of selectivity. 


Keywords: biosensors; electrochemistry; cellulose based membranes 





Introduction 


In recent years, considerable interest has been focussed on biosensors in 
which an enzyme or whole-cell coated sensor produces a signal in response to 
a substrate [1]. Many such systems utilise amperometric detection, and cov- 
ering membrane layers are important in this respect. Thus a semi-permeable 
membrane is often required between the enzyme layer and the electrode sur- 
face to furnish the desired degree of selectivity for the chosen analyte in bio- 
logical fluids. Avoidance by this procedure of complicated separation steps 
makes feasible invasive or “dipstick” sensors, suitable for clinical use [2]. 

A commonly used membrane is a dense cellulose acetate layer. This for ex- 
ample is utilised in a commercial glucose analyser [3], where H.O, is the de- 
tected end product of a glucose oxidase catalysed reaction. In this situation, 
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the inner selective membrane needs to be permeable to hydrogen peroxide, 
while excluding electrochemically interfering species from blood, notably as- 
corbate, urate, certain amino acids and small peptides, such as glutathione [4]. 

The successful use of such hydrogen peroxide selective membranes [4] in 
the glucose sensor has stimulated their incorporation into other oxidase based, 
direct responding electrodes. In reverse osmosis experiments, phenol and cer- 
tain of its derivatives are known to permeate dense cellulose acetate mem- 
branes more readily than inorganic species [5-7]. This offers scope for bios- 
ensor construction in situations where there is a need to detect a phenolic 
compound. 

A number of phenolic compounds are of clinical interest. Paracetamol ide- 
ally should be directly measured without the intervention of an enzyme [8- 
10]. Serum levels of both acid and alkaline phosphatase enzymes can be mon- 
itored by the dephosphorylation of phenolic compounds [11]. The resulting 
unconjugated phenols being electroactive, detection with an electrode is there- 
fore feasible, but as yet, not achieved in serum without a prior separation step, 
due to problems of background interference. 

Alkaline phosphatase is commonly used as a label in enzyme immunoassays, 
and electrochemical detection can also be applied here [12,13]. A number of 
substrates are suitable for the enzyme, e.g. phosphates of phenol and p-ami- 
nophenol, but these differ considerably in the ease of electrochemical product 
detection [13,14]. 

The work described here combines the phenol permeability of cellulose ace- 
tate membranes indicated by reverse osmosis, with the selectivity properties 
of previous biosensors. The responses of a polarised platinum electrode ob- 
tained through such cellulose acetate (CA) and cellulose acetate butyrate (CB) 
membranes are investigated. Comparison is made to the responses of a bare 
electrode and of one covered by an unselective microporous polycarbonate 
membrane (PCM). The phenolic species of clinical interest are considered 
along with a common biological interferent, ascorbate. Hydrogen peroxide and 
various other phenolics are investigated in order to correlate the structure and 
the permeation characteristics through the membrane. 


Experimental 


Materials 

Figure 1 shows the structures of the species investigated. Measurements were 
performed in isotonic phosphate buffer at pH 7.3, containing 2.44 g NaH.PO,, 
7.5 g Na,sHPO,, 3 g NaCl and 0.6 g EDTA per litre. A pH of 7.3 was chosen as 
being appropriate for oxidase enzyme electrodes and paracetamol measure- 
ments in serum. It was also adopted as the standard assay pH throughout the 
studies comparing the other species of interest. Polycarbonate membranes with 
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Fig. 1. Structures of species studied. 


track-etched pores of 2 4m diameter were obtained from Nuclepore® (Pleas- 
anton, CA). 


Preparation of membranes 

Solutions of cellulose acetate (39.8% acetyl content) and cellulose acetate 
butyrate (17% butyryl) were prepared by dissolving 1.5 g polymer in 100 ml 
acetone. One ml of the solution was evenly spread over a glass plate (5 x5cm7’). 
The glass plate was manually rotated, approximately once every two seconds 
during a two minute evaporation period; more rapid rotation led to uneven 
membranes being formed. Membranes were left to dry at room temperature 
for one hour. Membrane batches were prepared and stored for later study. Av- 
erage thickness, as assessed by a micrometer screw gauge, and by scanning 
electron microscopy in a limited number of membranes, was ca. 10 um. 
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Measurement procedure 

A two electrode cell was used, based on an oxygen electrode (Rank Brothers, 
Bottisham Cambridge, UK) (Fig. 2), with a platinum anode and silver cath- 
ode. For bare electrode measurements, the O-ring and membrane were omitted. 
Membranes covering the sensor were held in place by the O-ring to provide a 










































































1. Anode 5. Stirrer bar 

2. Cathode 6. Cell baseplate 

3. Membrane 7. Sample chamber 
4. O-ring 8. Locking ring 


Fig. 2. Cross section of the Rank electrochemical cell. 
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Fig. 3. Responses obtained to 0.04 mM solutions, using a bare electrode. 
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seal so that all responses were solely due to diffusion through the membrane. 
The platinum anode was polarized at 650 mV against the silver cathode, which 
acted as a pseudo reference, using an “in-house” potentiostat (School of 
Chemistry Workshops, University of Newcastle-upon-Tyne). This allowed 
control of polarising voltages (+ /—1.2 V) and detection of current flow in the 
range 0.1nA-2.0mA. A recorder output was connected to a Lloyd Instruments 
(Fareham, Hants) type CR652s strip chart recorder. 

The electrode was allowed to condition in buffer before use, until a stable 
base line was achieved (2 hr). Measurements were made by adding small vol- 
umes of stock buffered solutions to the assay buffer in the cell. A step change 
in current output resulted, and responses were recorded at the plateau value; 
the cell was rinsed with buffer between each measurement. 
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Fig. 4. Polycarbonate membrane; (a) responses obtained to 0.04 mM solutions, (b) responses 
using PCM, expressed as a percentage of the bare electrode responses. 














88 
Results 


The responses of the bare electrode to 0.04 mM solutions are shown in Fig. 
3. Both ascorbate and hydrogen peroxide were strongly electroactive, and se- 
lective membranes would be required to measure peroxide in samples contain- 
ing ascorbate. 

The relative responses are shown in Fig. 3. Phenol itself had the lowest elec- 
troactivity, compared to the other monohydroxy derivatives p-aminophenol 
and paracetamol. The dihydroxy compounds catechol, dopamine and CH1 also 
gave larger responses, while the trihydroxy species phloroglucinol was also 
electroactive, but less so than most catechol derivatives. 

When a PCM was used (Fig. 4a), the responses were reduced by approxi- 
mately 90% (Fig. 4b). Though the relative order of activities was unchanged, 
responses to the less active compounds were attenuated to a lesser degree. 

Cellulose membranes decreased the electrode signal to a greater extent. For 
both CA and CB membranes (Figs. 5a and 6a), peroxide responses were pro- 
portionally less reduced. The effects of cellulose membranes on the responses 
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Fig. 5. Cellulose acetate membrane; (a) responses obtained to 0.04 mM solutions, (b) responses 
using CA membranes expressed as a percentage of the bare electrode responses. 
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Fig. 6. Cellulose acetate butyrate membrane; (a) responses obtained to 0.04 mM solutions, (b) 
responses using CB membranes expressed as a percentage of the bare electrode responses. 


to different phenols were further assessed by determining the percentage of 
the residual signal with the membrane in place, compared to the bare electrode 
(Figs. 5a and 6b). This showed relatively little effect on the phenol signal, 
while those of the other species were greatly reduced. For a CA membrane 
phloroglucinol, catechol and particularly peroxide were less reduced than the 
other species. The results on CB were similar, but with all responses reduced 
to around half the value obtained with CA. The phloroglucinol, peroxide and 
catechol response was closer to that of the other compounds with this mem- 
brane. Figures 7 and 8 are calibration curves demonstrating a general linear 
dependence on concentration. This is critical for electrode calibration; the re- 
sults confirm that relative solute permeation through the membranes was con- 
centration independent, and up to 0.1 mM did not show saturation behaviour. 
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Fig. 7. Calibration curves obtained using cellulose acetate membranes. 
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Discussion 


Initial results were consistent with the observations of Tang et al. [13], and 
Koile and Johnson [15], that phenolic compounds can foul the surface of the 
electrode during oxidation. However, by adopting a low concentration (0.04 
mM), this problem was eliminated; no loss of signal in this study was experi- 
enced at this concentration. 

The electroactivity of the phenolics shown in Fig. 3 appears to show a rela- 
tionship to structure. The substituents -OH and —NH, have lone pair electrons 
adjacent to the ring, tending towards increased electron density in ortho and 
para positions, with the result that phenols with -OH or -NH, groups in ortho 
or para positions are more readily oxidisable. Hence the dihydroxy molecules 
catechol and dopamine and also the amines paracetamol and p-aminophenol 
gave larger responses than phenol. In the case of CH1, which is a substantially 
larger molecule, the relatively low response was due to its larger size which 
probably resulted in more restricted diffusion to the electrode surface. Phlo- 
roglucinol has three -OH groups, but these are in meta positions, which are 
less enhancing of functional group activity than either ortho or para 
configurations. 

The use of any membrane restricts diffusion to an electrode; the steady state 
concentration gradient at the electrode surface is correspondingly reduced, giv- 
ing proportionately lowered currents. 

A PCM brought about a general reduction in signal for all species. Such a 
membrane is microporous, and the effect can be attributed simply to a reduc- 
tion in area available for diffusion. 

The situation with cellulose membranes is somewhat different. These mem- 
branes constitute a homogeneous phase, and any permeation involves intimate 
solute contact with the membrane polymer. The distinction between homo- 
geneous dense membranes formed by a casting procedure similar to ours, and 
the Loeb-Sourirajan type with a dense surface layer, but porous bulk, has been 
demonstrated by electron microscopy [16,17]. The nature of the interaction 
will determine the flux through the membrane [5-7]. There are no large pores 
by which the diffusing molecules can avoid the influence of the membrane 
material. Cellulose membranes brought about a great reduction in signal size 
for all the species investigated. Hydrogen peroxide gave the greatest signal 
through both membranes and the ascorbate responses were amongst the low- 
est. This is ideal for use in oxidase enzyme electrodes, and the cellulose acetate 
membrane seems particularly suited. Reverse osmosis membranes of cellulose 
acetate were used by Matsuura and Sourirajan [18-20] with a range of phe- 
nolic and other organic species in aqueous solution, and a correlation of per- 
meation with pK, and hydrogen bonding ability was found. Nonionised mole- 
cules with a proton available for hydrogen bonding were more permeable than 
cations. Matsuura and Sourirajan [20] also described cellulose acetate as hav- 
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ing a proton acceptor character. Their results are in agreement with our data 
showing ascorbate to be excluded relative to the phenolics and hydrogen 
peroxide. 

A further general effect of cellulosic membranes is to reduce the high vari- 
ability of responses to various phenolics (Fig. 3), although catechol responses 
were the largest, and those to CH1, the smallest with both membranes. In 
particular, the already small phenol responses of the bare electrode were little 
reduced by the cellulosic membranes. For phenol, the low reactivity of the mol- 
ecule itself led to relatively small signals at the bare electrode, whereas for the 
membrane covered electrode, the decreased flux was a less drastic limiting fac- 
tor on signal size than for the more electroactive species. There was relatively 
little difference between the phenol derivatives. Charged dopamine, p-amino- 
phenol and CH1, and uncharged catechol, phloroglucinol and paracetamol did 
not fall into two readily discernable groups, the variation within each class 
being quite considerable. Cellulose acetate reverse osmosis membranes were 
found by Matsuura and Sourirajan to reject certain charged species, while al- 
lowing uncharged molecules to permeate [18,19]. They found that the cellu- 
lose acetate had a certain degree of proton acceptor character, which was re- 
sponsible for the rejection of anions. Since p-aminophenol and dopamine are 
cationic, their charge does not restrict membrane entry. In view of the anionic 
nature of CH1, the mechanism of its permeation is less clear, but is possibly 
explained by a partitioning effect, in which the anionic nature of the acid groups 
is offset by the uncharged part of the molecule. In addition to the rejection of 
ascorbate on the basis of charge, its non-aromatic nature may have been im- 
portant in diminishing membrane partitioning. 


Conclusion 


The behaviour of phenolic species on a bare electrode is directly related to 
their structure. Microporous membranes can act as a diffusion barrier, but 
without any selectivity properties. The use of cellulose based membranes can 
impact selectivity. 

For practical biosensors, it is the signal/interference ratio that is important. 
CA membranes are thus likely to be more appropriate than PCM or CB for 
selectivity, particularly for peroxide, against ascorbate (plasma _ [ascor- 
bate | <0.1 mM), which is the desired situation for oxidase based sensors. 

Similar results apply to p-aminophenol, phenol and paracetamol detection 
in the presence of biological fluids. It should however be noted, that parace- 
tamol remains a serious interferent in glucose determination with cellulosic 
membranes, and alternative polymers need to be developed for this problem. 
To exclude paracetamol by rendering it charged, would necessitate highly al- 
kaline conditions sufficient to hydrolyse cellulosic membranes. 
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Abstract 


The feasibility of using a membrane reactor to shift the equilibrium of formaldehyde production 
from methanol is investigated. Formaldehyde production from methanol over a silver catalyst is 
an endothermic, reversible, and conversion-limited dehydrogenation reaction. In the membrane 
reactor, the equilibrium is further shifted by removing the hydrogen through the membrane. This 
enables the membrane reactor to be operated economically at a lower temperature. The permea- 
bilities are measured using a porous Vycor glass membrane. The reaction kinetics are investigated 
over the temperature range of 300-400°C at atmospheric pressure in a fixed-bed differential re- 
actor over 99.998% silver needle catalyst. The membrane reactor performance is experimentally 
analyzed at atmospheric pressure to study the effect of changing operating variables such as tem- 
perature, space velocity, air/methanol feed ratio, and membrane surface area. A theoretical model 
is derived and a computer simulation is carried out based on the Runge-Kutta method to verify 
the above experimental data. Using an isothermal cocurrent model, a parametric study is con- 
ducted to provide basic information regarding scale-up and optimum operating conditions. The 
experimental study and parametric analysis of the present study point to the necessity of devel- 
oping an improved permselective inorganic membrane material. Such a new membrane will make 
the membrane reactor method an attractive alternative in many industrial applications. 


Keywords: glass membranes; membrane reactor; Vycor; formaldehyde production; porous 
membranes 





Introduction 


Equilibrium-limited conversion is common in chemical reactions. A number 
of industrial reactions show incomplete conversions because of equilibrium 
limitations. As these limitations seriously hinder or prevent commercial utili- 
zation of potentially useful reactions, there is enormous economic incentive to 
find technologies that can drive these reactions toward completion. A reaction 
is most simply driven to completion by removing the reaction products from 
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the reaction zone. Formaldehyde production from methanol over a silver cat- 
alyst is an endothermic, reversible, and conversion-limited dehydrogenation 
reaction [1-3]. In conventional formaldehyde production, the hydrogen is ox- 
idized with air to drive the equilibrium position resulting in increased produc- 
tion of formaldehyde. The conversion of methanol depends markedly on the 
temperature of the reaction [4-6]. The conversion improves as the tempera- 
ture increases, with the side-reactions (formation of carbon dioxide and car- 
bon monoxide) increasing at the same time. In the membrane reactor method, 
the equilibrium is shifted by removing hydrogen through the membrane. This 
enables the membrane reactor to be operated at a lower temperature, since the 
equilibrium shift induced by the removal of hydrogen through the membrane 
would eventually drive the reaction toward completion. At lower reaction tem- 
peratures, selectivity increases owing to the lesser extent of side-reactions and 
thermal decomposition of the reactant and the product. Therefore, the ex- 
pected low operating temperatures in the membrane reactor might provide a 
higher yield. 

The partial oxidation of methanol is the principal commercial route to form- 
aldehyde production. In general, this process is carried out by passing a mix- 
ture of methanol and air over a heated stationary catalyst at approximately 
atmospheric pressure and scrubbing the off-gases with water to obtain aqueous 
formaldehyde. The formation of formaldehyde from methanol and air was first 


believed to be a gas phase oxidation process, as indicated by the following 
equation: 


CH,OH+ 1/2 0.+CH,0+H,0 + 38kcal (1) 


This is still believed to be the mechanism followed when an oxide catalyst is 
employed. However, studies of the reaction mechanism and examination of 
the products obtained by the metal-catalyzed process demonstrated that under 
these conditions the mechanism is either dehydrogenation followed by oxida- 
tion of hydrogen to the extent that oxygen is present in the gaseous mixture as 


shown below, or a combination of the dehydrogenation and oxidation reactions 
[1,3,7]. 


CH,OH#@CH,0+ H, —20 kcal (2) 
H, +1/2 0. +>H,.O+58 kcal (3) 


Oxidation reactions supply heat to make the process self-sustaining and are 
also believed to keep the catalyst active and displace the dehydrogenation equi- 
librium to the right. Newton and Dodge [2] quantitatively determined the 
equilibrium relation for this reaction and found that 


K, = Pcu.0 Pu./Pcxson 


(4) 
=log-! (—4600/T+6.470) 
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The above experimental values of the equilibrium constants are compared with 
their thermodynamic values [8-10] and used in the calculation of the equilib- 
rium conversions. 

In recent years, many researchers have studied the feasibility of producing 
hydrogen by enhancing the conversion of dehydrogenation and decomposition 
reactions. Dokiya et al. [11], Fukuda et al. [12], and Kameyama et al. [13] 
employed a porous Vycor glass membrane to selectively separate the hydrogen 
produced in the decomposition of hydrogen sulfide. In a separate study using 
the same reaction, Kameyamaet al. [14] showed the possibility of using porous 
alumina membranes. Shindo et al. [15] conducted similar studies for the hy- 
drogen iodide reaction using a porous Vycor glass membrane. Ito et al. [16], 
in a theoretical study of the decomposition of hydrogen iodide, showed that a 
significant equilibrium shift is possible when a porous Vycor glass membrane 
is used to separate the hydrogen. Shinji et al. [17] and Ito [18], using a porous 
Vycor glass membrane and palladium, respectively, enhanced the dehydroge- 
nation of cyclohexane. Ito et al. [19] simulated the membrane reactor with the 
above reaction and showed that for given rates of permeation and reaction, 
there is an optimum thickness of membrane at which maximum conversion is 
obtained. In a separate parametric study by Mohan and Govind [20], it is 
known that for a membrane that is permeable to both products and reactants, 
the maximum equilibrium shift possible is limited by the loss of reactant from 
the reaction zone. Sun [21] studied the same reaction using a platinum-im- 
pregnated Vycor glass catalytic membrane reactor. A propylene disproportion- 
ation reaction has been studied by Chang [22] in a porous Vycor glass mem- 
brane reactor. This reaction is particularly interesting because it is reversible 
at room temperatures. Recently, Chen [23] studied the direct oxidation of 
ethylene to acetaldehyde in an aqueous solution of palladium (II) chloride- 
cupric chloride catalyst with a silicone rubber membrane reactor and a poly- 
propylene membrane reactor. 


Experiments 


Permeabilities of gases and vapors were measured using a porous Vycor glass 
permeation cell which was the same as the membrane reactor but without the 
catalyst and outlet. For gases (hydrogen, oxygen, nitrogen, and carbon diox- 
ide ), the pure gas was used as feed and the volumetric flow rate of the permeate 
gas was measured with a bubble flow meter. To measure the permeabilities of 
the vapor components, pure methanol, pure water, and a 37% formaldehyde 
solution (mixture of formaldehyde, water, and methanol) was vaporized in a 
separate heating vessel and each vapor component or a mixture of vapors was 
collected in a condenser using liquid nitrogen and a vacuum pump. 

The reaction rate constants were independently measured using a differen- 
tial fixed-bed reactor. The reactor was made small enough with a 1/4 inch 
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(O.D.) stainless steel tube to eliminate mass and heat transfer effects. The 
reactor was packed with 2.1 g of heat-treated 99.998% silver needle catalyst 
(mesh size 16x30). The catalyst was commercially available from Metz Me- 
tallurgical Corp., South Plainfield, NJ. The system operated at low conversion 
at atmospheric pressure so that the composition as well as the reaction rate 
could be considered constant at any point in the reactor. No appreciable cat- 
alyst deactivation was observed during the measurement of reaction rate con- 
stants between 300°C and 400°C. After initial operation for 1.5 hr to ensure 
attainment of the steady-state condition, the actual run was conducted for 1 
hr, during which the liquid and gaseous products were collected. The conver- 
sions were measured at various feed flow rates. 

A schematic diagram of the apparatus for the studies of the membrane re- 
actor is shown in Fig. 1. The membrane cell has a double-pipe configuration. 
The membrane is a porous Vycor glass tube (Corning Vycor brand glass, No. 
7930, code 742081, O.D. 10 mm, I.D. 8 mm) with one end connected to a 10 
mm ordinary Pyrex glass tube using a Vycor—Pyrex graded seal (Corning Cat. 
No. 6466). The other end of the inner tube is closed by glass fusion only. The 
porous section of the inner tube is ca. 8 to 21 cm in length. The outer shell is 
an ordinary Pyrex glass tube 16 mm outside diameter and 30 cm in length. One 
end of the inner tube is suspended in order to prevent breakage due to the 
uneven thermal expansions of the inner tube and the outer shell. No inlet line 
is installed in the permeate side chamber; thus the input molar flow to the 
permeate side is effectively kept at zero throughout the experiments. The tem- 
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Fig. 1. Schematic diagram of membrane reactor apparatus; (1) mass flow meter, (2) regulating 
valve, (3) metering valve, (4) syringe pump, (5) electric furnace, (6) porous Vycor glass mem- 
brane reactor packed with silver catalyst, (7) temperature indicator, (8) temperature indicator, 
(9) six-port valve, (10) GowMac 550 TCD gas chromatograph, (11) condenser with liquid nitro- 
gen, (12), ice-trap, (13) vacuum pump, (14) bubble flow meter, (15) pressure gauge, (16) plug 
valve. 
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perature is measured with three K-type thermocouples attached to the outside 
wall of the reactor cell, and is controlled within + 1°C of the desired value by 
an electric furnace and a Cole-Parmer digital temperature controller. The gas 
feed (dry air) is from a gas cylinder and the flow rate is measured with a mass 
flow meter. The liquid feed (methanol) is injected at a constant flow rate using 
a syringe pump, heated, and evaporated so that the vapor can be carried to the 
shell side (catalyst ) by the feed gas. Heating tapes are used to preheat the input 
stream and prevent condensation in the output stream. The outlet gas flow 
rate in the rejected stream is measured with a bubble flow meter, and the liquid 
outputs in both streams are collected in two separate condenser systems. A 
vacuum is applied to the permeate stream and feed side pressure is maintained 
at atmospheric pressure. A pressure gauge is used to measure the feed side 
pressure. The output stream was sampled with an on-line sampling valve and 
analyzed for the compositions of vapors (formaldehyde, water, and methanol ) 
using a Gow-Mac series 550 TCD gas chromatograph equipped with an Alltech 
6’ x 1/8” SS (80/100 mesh size) Porapak T column. Helium was used as the 
carrier gas (55 ml/min) and column temperature was 145°C. The composition 
of the exhaust gas stream (hydrogen, nitrogen, and oxygen) was determined 
using the same G.C. with an Alltech 3’ x 1/8” SS (65/100 mesh size) Molec- 
ular Sieve 5A column. Argon was employed as the carrier gas (25 ml/min) and 
the column temperature was 25°C. Carbon dioxide content was measured us- 
ing a Porapak T column at 25°C with a helium flow rate of 65 ml/min and 
compositions of carbon monoxide and methane were analyzed using a Molec- 
ular Sieve 5A column at 25°C with a helium flow rate of 35 ml/min. The feed 
composition of methanol and air is about 0.3 and 0.7, respectively, and the 
vapor compositions of the permeate and reject streams are checked every hour 
with the on-line sampling G.C. The liquid output flow rates are calculated 
based on the amounts of liquid collected after the experiment is finished. To 
obtain enough steady-state data, each experiment is allowed to last for 5 to 6 
hr. Just before the end of each run, the compositions and flow rates of the 
exhaust gases are to be measured. A plug valve is installed at the end of the 
permeate stream to allow comparison of the membrane reactor data with those 
of the packed-bed reactor. After each run, the valve is closed, and the packed- 
bed reactor performance is checked. 


Modelling 


For isothermal cocurrent model, the shell side of the membrane reactor is 
packed with a catalyst and becomes the reaction zone. The flow rate change at 
any intermediate point within the cell can be determined from a differential 
material balance as shown in Fig. 2. 
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Fig. 2. Cocurrent tube-and-shell membrane reactor. 


On the shell side: 

dF’ 

Gp ibe — 20Q Pua (xj —Pry;) An (R/Ri) (5) 
On the tube side: 

dG; 

Ge oA Pu (x) — Pry) In (Ro/ Ri) (6) 


The first term on the right-hand side of eqn. (5) represents the mass generated 
by reaction, where L, can be termed “catalyst linear density” with a measured 
value of 1.75 g/cm inside the tube. The second term accounts for the contri- 
bution from membrane permeation. Equation (6) implies that the flow rates 
on the tube side change solely due to permeation through the membrane wall. 
With the initial composition known, a fourth-order Runge-Kutta method is 
then used for the integration of the simultaneous differential eqns. (5) and 


(6). 


Results and discussion 


A. Permeability measurements 

The temperature dependences of the permeabilities in the porous Vycor glass 
membrane are plotted in Fig. 3 for the gases (hydrogen, oxygen, nitrogen, and 
carbon dioxide) and in Fig. 4 for the vapors (methanol, water, and formalde- 
hyde). The ratio of permeabilities of hydrogen and nitrogen in the studied 
temperature range (573—723°K) is 3.72, which is very close to the theoretical 
value of 3.74 predicted by the Knudsen flow mechanism in which the perme- 
ability ratio of two gases is inversely proportional to the ratio of the square 
roots of their molecular weights. Similarly, the permeability ratio of hydrogen 
to oxygen is 3.92, for which the theoretical value is 4.00. The permeability of 
carbon dioxide also conforms to the theoretical calculation with+3% error. 
However, the permeabilities of methanol, water, and formaldehyde do not fol- 
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Fig. 3. Gas permeabilities measured in porous Vycor glass membrane. 
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Fig. 4. Vapor permeabilities measured in porous Vycor glass membrane. 


low the prediction given by the Knudsen flow mechanism. For example, the 
experimental data for methanol deviate up to 39% from the predicted values 
at 573°K. The above deviations for those three vapor components could be 
explained by the mechanism of surface diffusion and condensation flow in this 
temperature range [24-26]. 
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B. Kinetic study 

Robb and Harriot [27] proposed the following kinetic equation, which is 
based on competitive adsorption of methanol and hydrogen on sites formed by 
reversible chemisorption of oxygen on silver and reaction of adsorbed metha- 
nol with underlying oxygen atoms: 


Pm = {RPp/(A+KmP mt KnPp) }{KoPe?/ (1+ KP?) } (7) 


This rate expression allows for a variable fractional order with respect to meth- 
anol and takes into account the effect of product inhibition and oxygen inlet 
partial pressure. Hydrogen is assumed to be the inhibiting product. If the ox- 
ygen adsorption constant is large, then the rate expression simplifies to 


Prm={RP,,/(1+KmPm+AKnPr)} (8) 


In order to obtain the constants in the rate expression the integrated form of 
eqn. (8) is fitted to the kinetic data using a nonlinear regression program 
(Marquardt’s ). 


kt’ =(K,P,3+1)In(1/(1—X) ) + (Kn — Ky) PmiX (9) 


A separate set of experimental data was gathered to see the effect of oxygen 
inlet partial pressure [28]. When P,,; is kept constant (0.95 atm) and the 
conversion level is very low, eqn. (7) can be reduced to 


rn=A {K,P9?/(1+K,Po)} (10) 


where A is assumed to have a constant value. The linearized form of eqn. (10) 
can be rewritten as: 


1/r,=1/A {1+ (1/K,) (1/P8?)} (11) 


From the plot of 1/r,, vs. 1/P°:°, K, is obtained. The results are summarized 
in Table 1. The values of k, K,,,, and K,, have the same order of magnitude with 
the literature value [27]. For the calculation of the activation energy, log k has 
been plotted against 1/7; the activation energy was found to be 10,942 kcal/ 
kmol. 


TABLE 1 


Values of rate constants 








Temp. kx 10° K,, K,, K, 

(°C) (kmol/kg-atm-sec ) (1/atm) (1/atm) (1/atm®*) 
300 0.2868 1.3 16.5 18.75 
350 0.5974 2.6 31.8 70.68 


400 1.1982 3.9 48.4 160.13 
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C. Performance of membrane reactor 
The membrane reactor performance was experimentally analyzed at atmos- 

pheric pressure to study the effect of changing operating variables. Atmos- 
pheric pressure was employed to maintain a high conversion level (fraction of 
reactant converted). Temperatures between 300°C and 400°C were chosen to 
maximize the extent of equilibrium shift and also to minimize any possible 
deactivation of the catalyst. Four different sizes of syringes were used to vary 
the space velocity (number of reactor volumes of feed treated in unit time). 

3 ml syringe; 0.009 ml/min, 0.3710 x 10~° kmol/sec 
10 ml syringe; 0.025 ml/min, 1.0306 x 10~° kmol/sec 
30 ml syringe; 0.057 ml/min, 2.3498 x 10~° kmol/sec 
50 ml syringe; 0.086 ml/min, 3.5453 x 10~° kmol/sec 
Conversion improves as the temperature and the space time (time required to 
process one reactor volume of feed) increase. In general, beyond a certain tem- 
perature the methanol conversion rises only slightly, whereas the side reac- 
tions show a more marked increase. The maximum yield (fraction of reactant 
converted into desired product) lies at the temperature at which the sum of 
the losses through the unreacted methanol and from the formation of carbon 
dioxide and carbon monoxide is at a minimum. Up to 400°C carbon monoxide 
is not detected and the formaldehyde selectivity increases. The yield increases 
with the increase in space velocity. This remarkable effect of the space velocity 
is apparently connected with the higher gas velocity in the catalyst bed. This 
may be explained by the fact that the turbulence of the gas molecules increases 
as they flow past the silver crystals, which increases the probability of their 
undergoing reaction [5]. The oxygen in the incoming air is almost completely 
consumed and a high hydrogen content must be maintained to obtain good 
yields at constant feed flow rate. This means that the dehydrogenation reac- 
tion should be encouraged as much as is consistent with practical conversion. 
When the temperature goes up and the space velocity increases, the water con- 
tent decreases. In the membrane reactor, the thermodynamic equilibrium is 
displaced by burning the hydrogen with air and further shifted by removing 
the hydrogen through the membrane. As expected, a longer space time gives a 
larger equilibrium shift as shown in Figs. 5, 6, and 13. 


D. Effect of air/methanol feed ratio 

The experiment was conducted at 60%, 70%, and 80% of air in the feed 
mixture to observe the influence of air/methanol feed ratio. The efficiency of 
the process depends to a large extent on maintaining an optimum ratio of 
methanol to air. Thomas [29] demonstrated that with the silver catalyst, gas- 
eous mixtures containing 0.5-0.3 gram oxygen per gram of methanol give bet- 
ter yields. A methanol—air mixture of the type cited contains 30-41 vol.% of 
methanol. Figures 5 and 6 show experimental results using a 12 cm long, 10 
mm Q.D. porous Vycor glass membrane reactor at a temperature of 350°C. 
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Fig. 5. Effect of space time on conversion at 60% air (350°C, 10 mm O.D., length = 12 cm porous 
Vycor glass membrane). 




















80 
Membrane Reactor (Simulation) 
70 4 . 
° 
Zz : Packed Bed 
ae 60 4 
¢ 
£ 
n 4 
— 
@® 
> 
e eed 
8 Equilibrium Conversion 
= 
40-4 
30 . + + T 
0.0 1.0 2.0 


Space Time (sec) 


Fig. 6. Effect of space time on conversion at 80% air (350°C, 10 mm O.D., length= 12 cm porous 
Vycor glass membrane). 


Comparison of these data with the computer simulation and the packed-bed 
reactor data is also included. At a higher percentage of air, the conversion and 
the selectivity increased while the equilibrium shift decreased. In other words, 
with increasing concentrations of methanol the manufacturing yields gradu- 
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ally fell off, but the net yield increased. Although the experimental data are 
scattered around the prediction given by the theoretical model depending on 
the space time, the isothermal cocurrent model simulates quite well the general 
trend of the actual membrane reactor system. 


E. Effect of membrane surface area 

To determine the effect of permeation area, four different lengths of mem- 
brane reactors were made with 10 mm O.D. porous Vycor glass membrane. 
10mm O.D., 8cm long porous Vycor glass; 22.3 cm? 
10 mm O.D., 12 cm long porous Vycor glass; 33.4 cm 
10 mm O.D., 17 cm long porous Vycor glass; 47.3 cm 
10 mm O.D., 21 cm long porous Vycor glass; 58.4 cm 
Figures 7 and 8 illustrate experimental results using two different sizes of sy- 
ringes at a temperature of 350°C. Higher conversions were obtained using longer 
membrane reactors. In Table 2, at a 70:30 air/methanol feed ratio 33.7%, 25.9%, 
and 16.0% of the maximum possible equilibrium shift was achieved at 300°C, 
350°C, and 400°C, respectively by the combination of the longest space time 
(6.0 sec) and the largest permeation area (58.4 cm”). In the above mentioned 
figures, the enhanced equilibrium shift by the membrane was 3.1%, 4.0%, and 
5.7% at 300°C, 350°C, and 400°C, respectively. It is apparent that the effect 
of the membrane reactor becomes greater at higher reaction temperatures as 
more hydrogen can permeate through the membrane because of its higher par- 
tial pressure. 

In this research, the catalyst linear density was fixed at 1.75 g/cm for all 
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Fig. 7. Effect of membrane surface area on conversion (350°C, 10 ml syringe, catalyst den- 
sity =0.63 g/cm’, 70% air). 
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Fig. 8. Effect of membrane surface area on conversion (350°C, 30 ml syringe, catalyst den- 
sity =0.63 g/cm’, 70% air). 


TABLE 2 


Equilibrium shift by membrane reactor (21 cm Vycor, 70% air, 3 ml syringe ) 








Temp. ( °C) , Xpb Xir 

300 23.1% 53.7% 56.8% 
350 43.3% 65.2% 69.2% 
400 65.8% 76.1% 81.8% 





membrane reactors. Due to this fact, the longer membrane reactor produced 
an adverse effect on the selectivity owing to the longer space time. For maxi- 
mum yields under this condition, a membrane length (or area) should be de- 
cided through the compromise between conversion and selectivity. Similar to 
the previous section, at longer space time (3 and 10 ml syringe) the experi- 
mental data were slightly lower than the model values. On the other hand, the 
shorter ones (30 and 50 ml syringe) had slightly higher conversions. This phe- 
nomenon implies that even if the outside wall temperature was controlled iso- 
thermally, the actual catalyst temperature at the larger space velocity (or 
shorter space time) might be somewhat higher than that of the smaller space 
velocity (or longer space time) due to the difference in the amount of air oxi- 
dation. For better control of the catalyst temperature, several researchers rec- 
ommended dilution of the feed mixture with water vapor [5,29-32]. The effect 
of the water can be explained as follows: some of the reaction heat from the 
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adiabatic process is absorbed by the water; in this way the oxygen fraction of 
the reaction mixture may be increased by passing in more air without the cat- 
alyst temperature rising. The dehydrogenation equilibrium is thus shifted in 
the direction of the formation of formaldehyde as a result of greater quantities 
of hydrogen being burnt. The yield rises to an optimum value as the proportion 
of water rises. If more water is added the yield drops in spite of increased con- 
version because more carbon dioxide is formed. 


F. Effect of membrane permeability 

The effect of membrane permeability was studied using the isothermal co- 
current model. Based on the experimental data of 12 cm long, 10 mm O.D. 
porous Vycor glass membrane reactor (catalyst linear density 1.75 g/cm), the 
permeation flux was increased up to 10 times. Results are shown in Fig. 9 for 
10 ml syringe and in Fig. 10 for 30 ml syringe at a temperature of 350°C. As 
can be seen in the previous section, increasing the membrane permeability (or 
the permeation flux) has the same effect as that shown by giving more mem- 
brane surface area. It is also observed that at higher temperatures the mem- 
brane reactor plays a more important role because of the hydrogen partial pres- 
sure. If the permeation flux is to be increased infinitely, the conversion reaches 
a maximum value and then drops gradually to zero since methanol partial pres- 
sure almost approaches zero value at some point due to the reaction and the 
permeation. In other words, the total conversion is attained by the competition 
of the hydrogen removal and the methanol depletion. Such information should 
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Fig. 9. Effect of membrane permeability on conversion (350°C, 10 ml syringe, membrane area = 33.4 
cm’, 70% air). 
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Fig. 10. Effect of membrane permeability on conversion (350°C, 30 ml syringe, membrane 
area = 33.4 cm”, 70% air). 


be taken into account in determining the permeation flux for a large scale 
operation. 


G. Effect of membrane selectivity 

The effect of membrane selectivity (separation factor or index of the relative 
separability ) was examined using the theoretical model. According to the ex- 
perimental data of 12 cm long, 10 mm O.D. porous Vycor glass membrane 
reactor (catalyst linear density 1.75 g/cm), the ideal separation factor, a, for 
this system is defined as 


where Q,, is the permeability of hydrogen through Vycor glass and Q,, is the 
permeability of methanol through Vycor glass. Results are shown in Fig. 11 for 
a 3 ml syringe and in Fig. 12 for a 10 ml syringe at a temperature of 350°C. 
Referring to the results in the previous section, it is found that increasing the 
membrane selectivity (or the ideal separation factor) gives the same effect as 
having more permeation flux. 


H. Recycle membrane reactor 

In a recycle membrane reactor a portion of the product is returned to the 
feed. As mentioned earlier, methanol is also permeable to the porous Vycor 
glass membrane and this will diminish the driving force for the forward reac- 
tion. If the permeated methanol is recycled, the above limitation of the mem- 
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Fig. 11. Effect of membrane selectivity on conversion (350°C, 3 ml syringe, membrane area = 33.4 
cm’, 70% air). 
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Fig. 12. Effect of membrane selectivity on conversion (350°C, 10 ml syringe, membrane area = 33.4 
cm’, 70% air). 


brane reactor system can be overcome. Figure 13 shows the comparison among 
the recycle membrane reactor (RMR), membrane reactor (MR) and packed- 
bed reactor (PR), at 350°C. All the values are taken from computer simula- 
tions using a 12 cm long, 10 mm O.D. porous Vycor glass membrane reactor 
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Fig. 13. Recycle membrane reactor (350°C, 10 mm O.D., length =12 cm Vycor glass membrane, 
70% air ). 


(catalyst linear density 1.75 g/cm). With the existing apparatus and operating 
procedures, it is difficult to recycle the permeated methanol, since the quantity 
is small and collected in a frozen state under vacuum and liquid nitrogen con- 
ditions. However, those difficulties may be resolved in a large scale process 
because commercial formaldehyde is obtained as aqueous formaldehyde by 
scrubbing the off-gas with water. As long as the recycle membrane reactor can 
be implemented into the real system, it will make up for the shortcoming of 
the membrane reactor method by a slight process modification. 


Conclusions 


The membrane reactor is experimentally analyzed and found to outperform 
the conventional process (the packed-bed reactor). The current study has not 
only demonstrated the enhanced equilibrium shift by removing the hydrogen 
through the membrane, but also introduced a new process using a membrane 
reactor as an economical alternative. Higher conversions are obtained using a 
longer space time and a larger membrane area. The effect of the membrane 
reactor becomes greater at higher reaction temperatures since more hydrogen 
can permeate through the membrane because of its higher partial pressure. 
The isothermal cocurrent model has been proposed and a computer simulation 
is made using the Runge-Kutta method to verify the experimental data. Al- 
though the experimental data are scattered around the prediction given by the 
theoretical model due to the heat effect, this model simulates quite well the 
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general trend of the actual membrane reactor system. The experimental results 
and parametric analysis of the present study point to the necessity of devel- 
oping an improved permselective inorganic membrane material, which should 
also be heat-resistant, inexpensive, and easy to handle. Such a new membrane 
will make the membrane reactor method an attractive alternative in many 
industrial application. 
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Subscripts 


eq 
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tube-side molar flow rate (kmol/sec ) 
shell-side molar flow rate (kmol/sec) 

forward reaction rate constant (kmol/kg-sec-atm ) 
adsorption equilibrium constant (1/atm) 
equilibrium constant (atm) 

catalyst linear density (kg catalyst/m reactor ) 
high pressure (Pa) 

low pressure (Pa) 

pressure ratio (=P; /Py,) 

permeability (kmol/sec-m-Pa) 

inside radius of membrane (m) 

outside radius of membrane (m) 

reaction rate expression (kmol/kg cat-sec ) 
temperature (°C or °K) 

conversion 

shell-side molar fraction 

tube-side molar fraction 

axial position (m) 

ideal separation factor (=Q,/Q,,) 

weight time (kg-atm-sec/kmol) 


equilibrium state 
hydrogen 

inlet condition 
component 
methanol 
membrane reactor 
oxygen 
packed-bed reactor 
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Absorbent Materials: Preparation, 
Structure, and Characterization. 1. The 





(J.R. Gross). 2. Preparation and structure 
of polyacrylates (F.L. Buchholz). 3. 
Preparation and characterization of 
crosslinked hydrophilic networks 

(L. Brannon-Peppas). 4. The equilibrium 
swelling behavior of porous and 
non-porous hydrogels(L. Brannon-Peppas 
and N.A. Peppas). 5. Characterization of 
ionic water-absorbent polymers: 
Determination of ionic content and 
effective crosslink density (D.C. Harsh and 
S.H. Gehrke). 6. Swelling behavior of 
water-soluble cellulose derivatives 

(E. Doelker). 7. Structure and swelling 
behavior of poly(ethylene 
glycol)/poly(methacrylic acid) complexes 
(J. Klier and N.A. Peppas). 8. Transport in 
porous polymers (W.M. Saltzman). 
Section Il. Recent Experimental 
Studies. 9. Molecular thermodynamics of 
aqueous polymers and gels (H.H. Hooper, 
H.W. Blanch, and J.M. Prausnitz). 10. 
Dynamic and equilibrium swelling 
characteristics of hydrophilic copolymers 
with carboxylic functional groups 

(A.R. Khare and N.A. Peppas). 11. Kinetics 
of swelling of absorbent polymers 

(FL. Buchholz). 12. Effect of 
superabsorber-pulp interactions on the 
performance of absorbent structures 

(D.L. Visioli and K.R. Williams).13.A 
microstructural analysis of polymer 
networks formed from graft copolymers in 
mixed aqueous solvents (C.G. Varelas and 
C.A. Steiner). Subject Index. 
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edited by M. Misono, Department of Synthetic Chemistry, Faculty of 


Engineering, The University of Tokyo, Tokyo, Japan, Y. Moro-oka, Research 
Laboratory for Resources Utilization, Tokyo Institute of Technology, 
Yokohama, Japan, and §. Kimura, Department of Chemical Engineering, 
Faculty of Engineering, The University of Tokyo, Tokyo, Japan 


(Studies in Surface Science and Catalysis, 54) 


The production of useful materials and 
the removal of polluting substances are 
fundamental to chemical technology, and 
in this respect catalytic and separation 
processes play essential roles. In order 
to cope with increasing demands to find 
solutions for the shortage of natural 
resources and global environmental 
pollution, rapid and significant progress 
in the technology is required. 


This book results from the successful 
seminar on Selective Reactions and 
Separation, held at Oiso, Japan, in 
February 1988. The seminar was 
organised by ASPRONC (the 
Association for the Progress of New 
Chemistry) as the fourth in a series of 
seminars on Frontier Technology. 
ASPRONC was inaugurated in 1986 and 
its membership comprises major 
companies in the chemical industry and 
various other sectors interested in 
chemistry. The aim of this seminar was 
to explore the frontiers of catalytic and 
separation technology and to discuss the 
requirements for its future development. 
The many interesting lectures and active 
discussions which resulted stimulated 
the editors to prepare this book. Each 
lecturer has written a chapter which 
represents a significantly revised and 
extended version of his original lecture. 


The book will appeal to many readers 
and will undoubtedly help to make a 
positive contribution to the future 
development of chemical technology. 





Contents: |. Present and Future of Catalytic and 
Separation Technology. Future opportunities in industrial 
catalysis (J.F. Roth). Design of practical solid catalysts at 
atomic/molecular levels (M. Misono). Characteristic dynamic 
behaviour of adsorbed species under reaction conditions 

(K. Tamaru). Recent advances in catalytic oxidations relevant 
to oxygenases (N. Kitajima, Y. Moro-oka). New trends in 
membrane separation technologies (S. Kimura). ll. New 
Aspects of Heterogeneous Catalysis. Novel catalytic 
functions of zeolites in liquid-phase organic reactions 

(Y. Izumi, M. Onaka). Activation and functionalization of 
methane (K. Otsuka). Catalytic decomposition of nitrogen 
monoxide (M. lwamoto). Effect of additives on the selectivity 
of metal catalysts (K. Kunimori et al.). Selectivity changes by 
alloying. Principles and application (A.J. den Hartog, 

V. Ponec). Recent advances and future developments in 
zeolite catalysis (Y. Ono). Ill. Catalysis Combined with 
Laser, Electrochemical and Membrane Technology. 
Application of lasers to control of chemical reactions - partial 
oxidation of alkanes in the gas phase (H. Tominaga, 

Y. Oshima). Electrocatalytic reactions with electron carriers 
(S. Torii). Recent development in membrane reactors 

(D.R. Seok, S.T. Hwang). Simultaneous operation of reaction 
and separation by a membrane reactor (N. Itoh). Functional 
membranes for separations based on reactions (T. Uragami). 
IV. Asymmetric Synthesis and Separation. New 
carbonylations by means of transition metal catalysts 

(I. Ojima). New chiral Rh(l) and Ru(!l) complexes: highly 
efficient catalysts for homogeneous asymmetric 
hydrogenation (H. Takaya et al.). Isolation of isomers and 
optical resolution by a host-guest complexation method 

(F. Toda). Recent progress and industrial application of 
optical resolution by high-performance liquid chromatography 
(K. Saigo). Author Index. Subject Index. 
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Thermodynamics of Polymer Solutions 
Phase Equilibria and Critical Phenomena 


by K. Kamide, Asahi Chemical Industry Co., Ltd., Osaka, Japan 


(Polymer Science Library, 9) 


This is the first self-contained book on the 
thermodynamics and critical phenomena of 
polymer solutions, ranging from the rather 
elementary level to the advanced and 
up-to-date level. The books covers the 
rigorous theories of phase equilibrium, 
computer experiments based on these 
theories, as well as actual experiments, 
molecular fractionation and application to 
membrane and fiber production. An 
extensive list of references and literature 
data on the thermodynamic interaction 
E-parameter, critical point, fractionation and 
polymer blends is also provided. This book 
should prove invaluable for courses on 
polymer science, thermodynamics and 
polymer solutions at graduate, university 
and polytechnic level. 


Contents. 1. Introduction. 2. Quasi-binary 
solutions of polydisperse polymers and a 
single solvent. Gibbs’ free energy of mixing 
AG. Chemical potentials A uo and A pyi. 
Two-phase separation characteristics. 


Fractionation based on solubility differences. 


Spinodal curve and critical solution point. 
Cloud point curve. Flory temperature and 
Flory entropy parameter. Flory enthalpy 
parameter. Experimental determination of 
phase diagram on polymer solutions. 

3. Quasi-ternary solution consisting of 
multicomponent polymer dissolved in 


binary solvent mixture. Chemical potential. 


Two-phase separation characteristics. 
Spinodal curve and critical solution point. 
Flory solvent composition. Cosolvency. 


Fractionation based on solubility differences. 


4. Quasi-binary solution consisting of 
multicomponent polymer 1 and 





multicomponent polymer 2 system. 
Introduction. Theoretical background of the 
critical solution point and spinodal curve. 
Theoretical background of the cloud point 
curve. Effect of polymer characteristics on 
two-phase equilibrium of multicomponent 
polymer 1 and multicomponent polymer 2. 
Experimental determination of 
polymer-polymer interaction parameter €1 2. 
5. Quasi-ternary solution consisting of 
multicomponent polymer 1, 
multicomponent polymer 2 and single 
solvent system. Introduction. Spinodal 
curve and critical solution point. Effects of 
molecular weight distribution of polymers 
and three —€ parameters on the spinodal 
curve and critical solution point. 

6. Application of phase equilibria: 
Formation of porous polymeric 
membrane by phase separation method. 
Introduction. Nucleation and growth of 
particles. Growth of primary particles by 
amalgamation to give secondary particles. 
Formation of pores by contacting the 
secondary particles. Critical points of casting 
solution/nonsolvent system and pore shape 
of membranes cast from the solution. 
Experimental relations between the 
solvent-casting conditions and pore 
characteristics of membranes. Formation of 
multi-layered structure of membranes. Some 
of the performances of multi-layered 
membranes. 7. Concluding remarks. 
Glossary of symbols. Author index. 
Subject index 
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Computer Applications in Chemical Engineering 


Proceedings of the European Symposium on Computer Applications in Chemical Engineering, 


ComChem 90, The Hague, The Netherlands, May 7-9, 1990 (21st Symposium of the Working 
Party on the Use of Computers in Chemical Engineering of the European Federation of the 


Chemical Engineering (418th Event)) 


edited by H.Th. Bussemaker, /BM Nederland, Zoetermeer, The Netherlands and 
P.D. Iledema, DSM Research B.V., Geleen, The Netherlands 


(Process Technology Proceedings, 9) 


This book presents the latest trends in computing, 
computer graphics and computerized design tools. It 
also gives a state-of-the-art overview of modelling, 
process integration and process design. All papers 
describe new computer algorithms and/or 
techniques for the whole range of computers from 
the PC to the supercomputer. Unit operations are 
well covered, as well as a number of topics in 
reactor engineering and control engineering. 


Contents. Session 1: Process modelling in revamp 
and design. Process integration - a new challenge to 
modelling. A dynamic simulator for design and 
analysis of chemical processes. A second generation 
nonequilibrium model for computer simulation of 
multicomponent separation processes. Process 
simulation with undefined mixtures using 
structure-based model compounds. Process 
integration retrofit subject to pressure drop 
constraint. Numerical methods used in reactor 
design (design of distributor top with internals) 
Using process synthesis and process simulation 
together - anew process engineering environment. 
Topology traps in evolutionary strategies for heat 
exchanger network synthesis. Modelling dynamic 
behaviour of multiple-effect falling-film 
evaporators. The use of sensitivity analysis to 
validate laboratory developed models of chemical 
reaction systems. Computer simulation for alcoholic 
fermentation process based on a heterogeneous 
model. MERITS, a simulation tool for multiphase, 
chemical reactor engineering. Session 2: Design 
and engineering graphics. Graphical 
representation of multidimensional functions in 
process design and optimisation. Design of 
integrated total flowsheet: a combined approach. An 
expert system for the design, operation and 
diagnosis of solid-liquid separators. User interfaces 
for process plant design - the benefits of Windows 
based systems. The architecture of SEPEX a 
SEParation system EXpert. Computer-assisted 
selection of equipment for liquid-liquid extraction. 
An algorithm for generating process flow diagrams 
automatically. Session 3: Process logistics and 





control. Models, modelling and management. 
Rigorous dynamic modelling and control system 
synthesis for distillation processes: an integrated 
approach. Control system design for an exothermic 
batch chemical reactor. Adding intermediate storage 
to noncontinuous processes. Scheduling ethylene 
plant operations. Simulation and optimization of an 
ethylene production plant. A computer-aided 
decision support system for the optimal long-term 
planning in large-scale chemical and related process 
systems. Computer aided analysis and control of a 
distillation train. Optimal operation of ethylene 
plants. On-line conversion estimation as an example 
of model supported process control. Session 4: 
Education and training. Education and training. A 
course in interactive graphics for undergraduate 
chemical engineers. Object-oriented hypermedia as 
a teaching aid in chemical engineering education. A 
simulation based system for operator training. The 
application of an operator training system in 
education. High fidelity real time simulation and its 
application to olefins plant operator training. 
Session 5: Trends in computing. Trends in 
computing. Computer aided molecular design by 
group contribution method. Flexibility, reliability, 
and availability analysis of manufacturing 
processes: a unified approach. A parallel computing 
strategy for multicomponent separation 
calculations. Optimal synthesis of nonsharp 
separation sequences. A mathematical framework 
for representation of heuristic rules. Two-phase flow 
simulation with cellular automata. Session 6: 
Information and data management in the 
chemical industry. Managing the IT function. A 
database approach for the interactive control of 
technical processes. Use of realtime and relational 
databases for CIM implementation. The case for a 
new database structure for computer-aided process 
design tools. Poster session. 
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An Introduction to 
RHeCOLOZY | rviicivresowesstonmtor 


excellent but mathematically rigorous 
texts on Rheology and a few on 
application in specific industries, there has 
been up to now no overview text suitable 
for the newcomer to the field that provides 
a painless introduction and overview to 





H.A. Barnes, Unilever Research, Port 
Sunlight Laboratory, Wirral, 

J.F. Hutton, formerly of Shell Research Ltd., 
Ellesmere Port, and 

K. Walters, University College of Wales, 


the subject and at the same time leads the 
Aberystwyth, UK i , 
panyeny reader with appropriate references to the 
(Rheology Series, 3) more detailed texts and literature 
references. 


Rheology is, by common consent, a difficult 
subject and some of the theoretical compo- 
nents are often viewed as being of prohibitive 


“Barnes and Hutton, with their years of 
industrial experience in application of 


complexity by scientists without a strong Rheology, have combined with the 
mathematical background. There are also the academic maturity of Walters to create the 
difficulties inherent in any multidisciplinary first text for the beginner in Rheology - 
science like rheology for those with a specific either for the student who has not had the 
training. Therefore, newcomers to the field are appropriate formal training, for the 
sometimes discouraged, and for them the industrial scientist, or for the industrial 


existing texts on the subject are of limited 
assistance because of their depth of detail and 
highly mathematical nature. 


manager who needs sufficient information 
to be convinced that knowledge extending 

beyond a shear rate dependent viscosity 
This book introduces the subject of rheology in may be necessary to solve his or her 
terms understandable to non-experts. It takes problem... Overall, the book will be a 
the seadier step by siep through the subject, and necessary part of your Rheology library, 
describes the application of rheological b : : ; 

ut most importantly it provides a well 


inciples to many industrial products and ; ; : 
as It pen Po norm — _— written and useful introduction to the 


authoritative guide which shows clearly how subject and thus will be the first book for 
mathematics, physics and chemistry have the newcomer to the subject.” 
contributed to the development of rheology. David V. Boger, The University of 
The generic features of all liquid-like materials Melbourne, Australia 








are summarised, i.e. viscosity, linear 
viscoelasticity, normal stresses and extensional 
viscosity. Particular systems are discussed, i.e. Contents: 1. Introduction. 2. Viscosity. 3. Linear 
polymeric liquids and suspensions. The final _ Viscoelasticity. 4. Normal Stresses. 5. Extensional 





chapter gives an outline of the theoretical Viscosity. 6. Rheology of Polymeric Liquids. 7. 
advances which have been made. Consistent | Rheology of Suspensions. 8. Theoretical Rheology. 
notation and nomenclature have been used Glossary of Rheological Terms. References. Author 


throughout the book, and the key textbooks Index. Subject Index. 
and publications which will enable the reader 
to follow up particular topics are listed. 
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